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SUMMARY 


An empirical method for the determination of the area, rate, and 
distribution of water-drop impingement on airfoils of arbitrary section 
is presented. The procedure represents an initial step toward the 
development of a method which is generally applicable in the design 
of thermal ice— prevention equipment for airplane wing and tail sur- 
faces. Results given by the proposed empirical method are expected 
to be sufficiently accurate for the purpose of heated-wing design, 
and can be obtained from a few numerical computations once the velocity 
distribution over the airfoil has been determined. 

The empirical method presented for incompressible flow is based 
on results of extensive water-drop trajectory computations for five 
airfoil cases which consisted of 15— percent— thick airfoils encompassing 
a moderate lift -coefficient range. The differential equations per- 
taining to the paths of the drops were solved by a differential 
analyzer. 

The method developed for incompressible flow is extended to the 
calculation of area and rate of impingement on straight wings in sub- 
sonic compressible flow to indicate the probable effects of compressi- 
bility for airfoils at low subsonic Mach numbers. 


INTRODUCTION 


The design of thermal ice— prevention equipment for airplane wing 
and tail surfaces has progressed to the point where the amount and 
distribution of heat flow can be calculated for specified flight and 
icing conditions (reference l) . This design procedure requires infor- 
mation as to the area, rate, and distribution of water-drop impinge- 
ment on the leading edge of the airfoil section being analyzed. At the 
present time, an approxi m ation of area and rate of water-drop impinge- 
ment is achieved by using a method involving the substitution of a 
circular cylinder for the airfoil leading edge, as suggested in 
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references 1 and 2. This substitution method is adequate for design 
purposes for some combinations of cylinder diameter and drop size, 
but it can produce sizable errors for other combinations (references 
1, 3> and 4) . 

A second means of estimating the area and rate of water-drop 
impingement on airfoils is provided by reference 3. This method is 
more accurate than the cylinder substitution method, but the calcula- 
tion procedure is somewhat laborious and, as a result, its use is not 
too practicable in a complete design study where a large number of 
trajectories are usually required. 

To establish a procedure which would eliminate the laborious 
computations of water-drop trajectories in the design of wing thermal 
ice— prevention equipment, it became apparent that a large number of 
water-drop trajectories would be required for study. Experience with 
calculating trajectories by the method of reference 3 had shown that 
the pattern of water-drop impingement for drop sizes usually encoun- 
tered in flight can be related most directly to velocity distribution 
over the surface of the airfoil. Airfoil shape itself appeared to 
have an effect on the pattern of impingement, but to a lesser degree 
than velocity distribution. Five airfoil cases were chosen as being 
the minimum which could be expected to provide sufficient data to 
include the effects of these two factors. Water-drop trajectories 
were computed for these five cases. 

This- report presents the results of the water-drop— trajectory 
computations. From a generalization of these data, a method is 
derived that permits a relatively rapid determination of the area, 
rate, and distribution of water-drop impingement for a fairly large 
number of airfoil surfaces of arbitrary profile. The derivation 
and limitations of this method and the procedure for its use are 
presented herein. The water-drop trajectories for the five airfoil 
cases were computed on the differential analyzer at the University of 
California, Los Angeles, under contract with the NACA; and appreci- 
ation is extended to the staff of the University for its aid and 
cooperation during the conduct of the program. 


, / 

SYMBOLS 


The following nomenclature is used throughout this report: 

a airfoil mean— line designation, fraction of chord from 

leading edge over which design load is uniform 

a d instantaneous drop-acceleration ratio, dimensionless 
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area normal to flow direction outlined by several trajectories 
at free— stream conditions, square feet 

area of impingement outlined on an airfoil surface by 
trajectories starting at free— stream conditions from an 
initial area of A 0 , square feet 


chord length of airfoil, feet 

/dA 

concentration factor ( • 


o 

V M S 


)' 


dimensionless 


drag coefficient of drop, dimensionless 

section lift coefficient, dimensionless 

'U \2 


pressure coefficient 


[ 1 


dimensionless 


pressure coefficient at Mach number M, dimensionless 


collection efficiency of airfoil, percent 

rate of change of velocity along the stagnation stream- 


line at the stagnation point 


" d(U a /v) 

“ dS 



dimensionless 


frontal height of airfoil, fraction of chord 

slope of airfoil contour at a particular chordwise 
position, dimensionless 

length of span, feet 

liquid-water content of icing cloud, .pounds of water 
per cubic foot of air 


Mach number, dimensionless 

weight rate of water-drop impingement per unit of surface 
area, pounds per hour, square foot 

weight rate of impingement of water drops on a body, per 
unit span, pound per hour, foot 


P ratio of the vector difference between the local air and 

drop velocities to free— stream velocity 
dimensionless 
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t t 

x ,y 


Ay 


radius of drop, feet 
Reynolds number for drop at 


Reynolds number for drop at 

distance along airfoil surface from leading edge, positive 
on upper surface and negative on lower surface, feet 

distance along water-drop trajectory, fraction of chord 

time, seconds 

equivalent ellipse thickness ratio for a low-drag airfoil 

( t— — — ), fraction of chord 

\ T max / 

maximum thickness of airfoil, fraction of chord 

component of local velocity parallel to chord line, 
feet per second 

local velocity of air or drop, feet per second 

component of local velocity perpendicular to chord line, 
feet per second 

free— stream air velocity, feet per second 

rectangular coordinates for a system of axes having the 
origin at the airfoil leading edge and the x axis, 
positive toward the trailing edge, lying along the 
airfoil chord, fraction or percent of chord 

rectangular coordinates for a system of axes having the 
origin at the airfoil leading edge and the x* axis, 
positive in the free— stream direction, lying parallel 
to free— stream direction, fraction or percent of chord 

total airfoil— ordinate intercept established by two 
impinging trajectories starting from infinity at a 
distance Ay Q apart, fraction of chord 

distance between two trajectories at infinity, fraction 
of chord 


relative velocity F7 
free— stream velocity Y 


(¥) 
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distance between two trajectories which start at infinity 
and impinge tangentially on the airfoil, fraction of 
chord 


angle of attack, degrees 
Prandtl-Glauert factor, equal to 




1/2 


) 


dimens ionle s s 


specific weight, pounds per cubic foot 


angular displacement between local velocity and -x axis, 
degrees 

kinematic viscosity of air, square feet per second 
airfoil leading— edge radius, fraction of airfoil chord 


time scale 


scale modulus) 9 


(?) ' « 


dimensionless 


, dimens ionle s s 


stream function, dimensionless 


Subscripts 

air 

average 

critical 

drop 

effective 
lower surface 

conditions at a particular Mach number 
maximum 

initial condition 
condition at airfoil surface 
tangential 
upper surface 
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^DERIVATION OF THE METHOD 

The method, derived herein, for calculating area, rate, and 
distribution of drop impingement assumes that airfoil velocity distri- 
bution can be considered as the primary factor influencing the paths 
of water drops which approach an airfoil. Hence, it was desired to 
find a relatively simple relation between the airfoil velocity distri- 
bution (which in itself is a function of airfoil contour) and the area, 
rate, and distribution of impingement. Accordingly, trajectories 
around airfoils of known surface velocity distribution were obtained 
for analysis, and, from this information, generalizations were made 
as to the behavior of drops in the presence of the different airfoil 
velocity fields. 


Description of Procedure Used to 
Obtain Water-Drop Trajectories 


The five airfoil cases selected for the water-drop— trajectory 
investigation are listed in the following table: 


Case 

Airfoil 

Angle of attack 
(deg) 

c l 

Leading— edge 
radius 

(percent chord) 

1 

15— per c ent — t h i c k 
symmetrical 
Joukowski 

0 

0 x 

2.67 

2 

Do • 

2 

.22 

2.67 

3 

Do 

4 

.44 

2.67 

4 

15— per c ent—' th ic k 
cambered 
Joukowski 

0 

.44 

2.67 

5 

MCA 65 2 -015 

( symmetrical) 

4 

.44 

1.505 


This table shows the systematic changes in the variables which 
affect velocity distribution. Cases 1, 2, and 3 were intended to 
reveal the effects of systemmatically altering airfoil velocity 
distribution by changing angle of attack; case 4, compared to cases 1 
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and 3, the effects of altering velocity distribution by the addition 
of a basic load distribution obtained by cambering the mean line; 
and cases 3 and 5, the effects of changing general airfoil shape for 
a given angle of attack and lift coefficient. The upper— and lower- 
surface velocity distributions over the forward region of each of the 
five airfoils are shown in figure 1. Velocity distributions for 
Joukowski airfoils are used where possible because the required veloc- 
ity components in the field of flow are more readily obtained than for 
most other airfoils. It is noted in figure 1 that the variables 
selected did not result in a wide variety of velocity distributions, 
but it is believed that these distributions are of sufficient scope to 
allow a generalization of the trajectory data, at least for cases in 
which there are no marked nose— pressure peaks. 

. The water-drop— trajectory computations were made to encompass 
a speed range of 100 to 350 miles per hour (assuming incompressible 
flow), a drop— diameter range of 20 to 100 microns, and a variation 
in altitude from 'sea level to 20,000 feet. Airfoil chord length was 
varied from 3 inches to 30 feet. These variables were combined into 
the dimensionless parameters, and Ry, which then were used as 
the independent variables throughout the trajectory computations. 

The range in values of t and Ry resulting from a combination of 
each minimum value and a combination of each maximum value of the three 
constituent variables is about 150 to 20,000 for and about 35 to 
1000 for Ry. These ranges in 'l' and Ry encompass most other pos- 
sible combinations of the selected values of speed, drop size, altitude, 
and chord length, but the data calculated for these ranges in \|r and 
Ry are not necessarily limited to the particular values of the con- 
stituent variables used to establish the ranges in \|/ and Ry. 

The problem of obtaining area, rate, and distribution of water- 
drop impingement on an airfoil is one of determining the solution to 
a set of simultaneous differential equations yielding the trajectory 
or path which a water drop will follow. These equations, a derivation 
of which may be found in reference 5, are essentially those which result 
from imposing conditions of dynamic equilibrium on a drop moving in an 
air stream. In dimensionless form, the equations are 
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Basically, equations (l) and (2) define the acceleration of a 
drop at any instant in orthogonal (x and y) directions. Consequently, 
a double integration of these equations, starting from a selected 
initial point (x 0 > yo)> yields x and y values necessary to plot a 
drop trajectory. Equation (3) is required to proceed with the solu- 
tions of equations (l) and (2). In performing the integrations, 
knowledge of the quantity CaR/2'4 (the ratio of the' actual drag coef- 
ficient given by Stokes law of resistance) is required; also required 
are magnitudes of the air— velocity components u a /V and v a /V as a 
function of drop location relative to the body. (See reference 5.) 
Variation of the term with local Reynolds, number R was taken 

from reference 6, while the variation of the air— velocity components 
u a /v and v a /V throughout the flow field was obtained analytically 
for the Joukowski airfoils. In the case of the NACA 65 2 -015 airfoil, 
however, the velocity distribution throughout the flow field was 
obtained by an electrolytic analogy method. 1 

In carrying out the differential analyzer computations for the 
five airfoil cases, the general procedure was to assign values to the 
terms t and Rv in equations (l), (2), and (3), to establish initial 
conditions, and then to obtain the water-drop— trajectory traces from 
the analyzer. For each combination of i|r and Ry selected, several 
trajectories were traced until the two trajectories were fo un d, one 
for the upper surface and one for the lower surface, which were t an gent 
to the airfoil surface at the point of drop impact. The importance of 
these two tangential trajectories lies in the fact that all drops 
between the tangential trajectories hit the airfoil and all drops out- 
side will miss. In some cases, after the tangential trajectories were 

electrolytic analogy method was based on the fact that, the stream- 
lines in an inviscid incompressible fluid and the equipotential 
lines in an electrical field are governed by the same equations. By 
means of this analogy and suitably constructed apparatus, velocities 
at any point in the flow field around a body could be measured 
directly. Accuracy of the electrolytic analogy method of obtaining 
velocities was assessed by comparing measured, velocities with theo- 
retically calculated values for two bodies. Comparisons were made 
for a small cylinder, with and without circulation, and for a 15- 
percent-thick Joukowski airfoil (of various chord lengths) with and 
without circulation. For optimum test conditions, the accuracy in 
velocities was about ±5 percent out to approximately one-half chord 
forward of the airfoil. Between one-half and one chord length 
forward of the airfoil, velocities could be obtained with an esti- 
mated accuracy of ±10 percent. 
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established, the distance between them was divided into six approxi- 
mately equal spaces, and trajectories started at the boundary of each 
space were traced. These intermediate trajectories were used to 
obtain an indication of the distribution of water-drop impingement 
over the airfoil surface. 


Water-Drop-Trajectory Data 


All the tangential trajectories, and the trajectories traced for 
distribution of impingement, are presented for the five airfoil cases 
in figures 2 through 6. The. drop trajectories, which were calculated 
for assigned values of the independent variables \|r and By directly 
yield values of trajectory starting ordinates and surface positions 
of drop imp ing ement. These variables, which, by the method of analysis 
shown in reference 1, can be used to determine the values of the 
dependent variables, area, rate, and distribution of impingement, are 
tabulated in tables I through V for each of the airfoil cases. Also 
included in the tables are the velocity components of the drops at the 
point of impingement. 


Manner in Which Water-Drop— Trajectory Data 
Were Generalized 


The objective in generalizing the water— drop— trajectory data 
obtained from the differential analyzer is to establish a relatively 
simple method for determining the impingement pattern of water drops 
on an arbitrarily selected airfoil at a given attitude for specified 
values of ^ and By. It is desirable that application of such a 
method should require only information which is usually available for 
the airfoil profile, or which is easily procured, such as the pressure 
distribution. Therefore, the first part of this section will present 
the reasoning associated with deriving such a method for incompressible 
flow from the results of the trajectory data; the latter part is 
devoted to an extension of the basic method to subsonic compressible 
flow . 


Generalization in regard to area of water-drop impingement .— In 
order to determine the area of water-drop impingement on the leading 
edge of an airfoil for specified meteorological and flight conditions, 
the values of s/c for the trajectories which impinge tangentially on 
the upper and lower surfaces must be obtained. In computational 
methods like those of references 3, 5, and 6, the procedure essentially 
has been to select values of ' l l' and By and then to .determine the 
trajectory. Various trajectories are computed until the tangential 
trajectory for the upper and lower surfaces is found. The two 
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tangential trajectories located determine the farthest positions of 
drop impingement on the airfoil surface for the selected values of t 
and Ry and permit calculating area of impingement from the equation 



In the derivation of the empirical method, herein presented, the 
reverse procedure was employed; that is, a point on the airfoil was 
selected (s/c) and the corresponding ^ and Ry values which are 
associated with the tangential trajectories at that point were deter- 
mined. The nature of the relationship between s/c and the parame- 
ters \|r and Ry is shown in figure 7. Data for the figure are those 
of table IV for the cambered airfoil at zero angle of attack and a 
lift coefficient of 0.44. From figure 1 , it can be seen that any 
specified value of s/c in the figure can correspond to an infinite 
number of combinations of Ry and Consequently it becomes 

necessary to select values of one variable and to solve for the other. 
In the derivation of the method of this report, values of Ry are 
assumed and corresponding values of are computed. 

If the data of figure 7 could be made available for all airfoils 
of interest, the problem of determining s/c for various values of 
+• and Ry would, of course, be simple. Because obtaining such data 
for all airfoils is impractical, the problem in the general case 
arises in determining the relation between Ry, \Jr , and s/c . To 
determine this relationship, equations (l), (2), and (3) are utilized 
to derive an expression for the scale modulus V , which is 


t = 


a d E V 



(4) 


where 



Equation (4) expresses generally the relation between ^ and Ky 
at all points in a trajectory, and, therefore, it is applicable at the 
airfoil surface for an arbitrarily selected value of s/c which corre- 
sponds to some particular tangential trajectory. It remains to 
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establish the values of C^r/24, R/fiy, and. a^ for the selected value 
of s/c. Actually, since C<iR/2k is a known function of R, the 
problem reduces to approximating R/Ry and a<i at the airfoil sur- 
face. 

Evaluation of R/Ry at airfoil surface: To determine R/Ry 

the method of this report is based on a graphical solution utilizing 
the hodograph plane. It is in this connection that the trajectory 
data from the differential analyzer are plotted in the hodograph 
plane in figure 8 for each of the five airfoil cases. To show the 
general relation of drop velocities to air velocities, the hodograph 
of air at the airfoil surface is also shown for each case. It is 
interesting to note in figure 8 that the velocity components for all 
drops, regardless of the combination of \|r and Ry, can be repre- 
sented by one faired curve. 2 Of particular interest, however, is the 
fact that the hodograph for the drops, for both upper and lower air- 
foil surfaces, always passes through the point u^/V = cos a, 
v^/V = sin a. Thus, in the simplest case of an airfoil at zero 
angle of attack, the hodograph of the drops always passes through 
an abscissa value of unity because this point corresponds physically 
to the point of maximum airfoil thickness where straight-line tan- 
gential trajectories always impinge upon the airfoil with free— stream 
air velocity. The coordinates at the origin of the air and drop 
hodographs correspond physically, of course, to the airfoil stag- 
nation point. 

To show the connection between the results of the trajectory 
data for the five airfoil cases and the case of an arbitrary airfoil, 
figure 9 is presented. Figure 9(a) depicts several water-drop 
trajectories in the physical plane impinging tangentially at the 
same point s/c on an airfoil which is at an angle of attack a. 
These trajectories, for purposes of explanation, may be assumed to 
correspond to various points on a line of constant s/c in figure 7. 
Accordingly, there is an infinite number of particular combinations 

2 An inspection of figure 8(e) will show that, for the upper surface 
of the MCA 65 2 -015 airfoil, the curve was not faired through all 
the data points. The reason is that the velocity components for 
some of the data points were in error. This can be shown by noting 
that every data point will not fall on a line connecting the origin 
nnfl a point on the air hodograph having the same s/ c value as the 
data point. Each data point must lie on such a line because the 
slope of a trajectory tangential to the airfoil at some s/ c posi- 
tion nnri the slope of the airfoil at the same s/c position are 
equal. The data points in figure 8(e) considered to be inconsistent 
with this criterion are designated with tag marks" and are not 
heavily relied upon for establishing the drop hodograph. 
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of ^ and Ey which are affine to any particular position of tan- 
gential drop impingement (s/c)^. In figure 9(a), a single vector 
representing the drop velocity for all the trajectories is drawn 
tangentially to the airfoil at the point of drop impingement . Only 
one vector is shown because the tangential trajectory hodographs 
presented in figure 8 indicate that all, drops impinging tangentially 
at a common point may be considered to have the same velocity. Also 
shown in figure 9(a) is a vector representing the air velocity at 
the point of impingement for the trajectories. The angle between 
the drop— and air— velocity vectors and the x axis is designated by 
the angle 9 . In figure 9(b), a typical air and drop hodograph is 
shown and the same vectors as shown in the physical plane are indi- 
cated. It should be noted that the difference in length of air and 
drop vectors at a particular s/c position is numerically equal to 
the value of E/By given by equation (3). 

As can be seen from figure 8, all five airfoils exhibit the same 
characteristic in regard to the single— valued nature of the term 
B/By with chordwise position, and it will be assumed that all other 
airfoils will display the same characteristic. However, in order to 
calculate values of B/By for an arbitrary airfoil, both hodographs 
of the air and of the tangential trajectories are required. The 
hodograph of the air velocity at the airfoil surface is easily 
obtained from the velocity distribution over the airfoil, but the 
shape of the hodograph for the tangential trajectories is more dif- 
ficult to obtain. This is because the data of figure 8 yield only 
the information that the tangential— trajectory hodograph always will 
pass through the point u^/V = 0, v^/v = 0 and the point 

= cos a, v^/v = sin a. 

With two points on the trajectory hodograph always known, it was 
concluded that, if one more point could be established, preferably 
where the vertical— velocity component reaches the maximum value, 
the general shape of the trajectory hodograph would be reasonably 
ascertainable. It was noted from the hodographs in figure 8 that 
peak values of v a /v and v^/v were at nearly the same location on 
the airfoil surface. That is, values of v a max /V and v d max /V' 

nearly always fall on a straight line through the origin, and it 

remains to determine just where on the line the value of vh /v 

max 

will fall. In this connection, it was decided to compare, for the 
five airfoil cases, values of the vertical component of relative 
velocity between drop and air attained at the position of maximum 
vertical air velocity. Accordingly, values of ( v a mflY A) — (v^ niay /V) and 

v amax /V were obtained from figure 8 and these are plotted in 

figure 10. An inspection of the data in figure 10 shows that the 
four Joukowski airfoil cases provide data sufficient to establish a 

relation between (v a /V) — (v^ /v) and v a /V; but too little 

max * max max 



NACA TN 2kj6 


13 


information is available for the NACA 65-series sections to firmly 
establish a relation between (v^A) ~( v d msT ^ and V a max / v * 

The point plotted in figure 10 for the NACA 65 2 -015 airfoil upper 
surface does not }.ie on the curve established by the Joukowski air- 
foil data, flnrl a question arises as to whether this difference is 
real or due to a possible improper fairing in figure 8(e) of the 
upper— surface drop hodograph. While this question cannot be resolved 
until further data are available, qualitatively it would seem that 
the drop velocities should tend to approach more nearly the air 
velocities in the case of low-drag airfoils because these shapes 
are not so conducive to abruptly altering the paths of water drops. 


To show how figure 10 can be used as an aid in the construction 
of the drop hodograph for an airfoil of arbitrary profile and 
attitude, figure 11 is presented. The procedure is first to draw 
the air hodograph and establish Vg^/v. Then, the maximum vertical 

velocity of the tangential— trajectory hodograph is determined as 
being less than V a max / v the amount ( v a fflax / V ) ~ ( v d Tnay / V ) in 

accordance with the curve in figure 10. The value of Y & rmx ^ so 


determined is assumed to lie on a straight line connecting the 

v a /V. The position of v d /V along the radial line 
max 


origin and 
determines the value of 

R/R^ 

s/c 


(R/%) 


‘V'v c 


^max 

at that particular position. 


a max 

Values of R/Ry for other s/c positions are taken, as a first 
approximation, as being in the same ratio to the air velocity at the 
particular s/c position as the value of R/Ry at v a /V is to 


U a /V at v^/V (curve A in figure ll). Values of R/Ry calcu- 
lated in this manner usually are too large near the point 
Ua/V = cos a, v d /v = sin a (point X, figure ll), so that a drop 

hodograph so constructed would probably not pass through this point 
as it should. To overcome this discrepancy in the drop hodograph 
as computed by a proportion based on the peak point of the air 
hodograph, a curve without reflex is faired tangentially into the 
proportional drop hodograph from the point u^/v = cos a, 
v d /V = sin a,. The combination of the proportional curve and the 

faired curve comprises the drop hodograph, which is labeled curve B 
in figure 11. For the five airfoil cases, maximum deviations of 
the drop hodographs obtained by the foregoing method from the drop 
hodographs in figure 8 were of the order of 15 percent. 


Two other methods were considered for establishing drop hodo— 
graphs. One method assumed R/Ry to maintain a constant value equal 

to the value prevailing at the point u d /v = cos a, v d /v = sin a. 


I 
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The other method assumed the ratio 



to maintain a constant 


value determined by the value of R/Ry and U a /v at the point 
Ud/V = cos a, v d /V = sin a. The drop hodographs given by each of 
these two methods also are shown for the example in figure 11. The 
curves are labeled C and D, respectively. These two methods have 
the advantage of not requiring the use of the hodograph and fig- 
ure 10; however, they are considerably more inaccurate (maximum 
deviations from the drop hodographs for the five airfoil cases being 
in the order of 30 percent) due to the neglect of factors of apparent 
influence on the drop trajectories. Use of either one of these 
latter two methods is suggested for particular airfoil cases which 
might happen to fall beyond the scope of the data presented in fig- 
ure 10. 


After the tangential— trajectory hodograph has been established 
in relation to the hodograph for air, values of R/Ry are available 
for various chordwise positions on the airfoil. These values are used 
in equation (4) for arbitrarily selected values of Ry and s/c. 

Since values of Ry are arbitrarily selected, values of R are 
ascertainable. Furthermore, values of the term / C d R/24 can be calcu- 
lated because Q^R/2k is a function of R in accordance with 

table VI. Thus, to solve equation (4), the only other term to be 
evaluated is a d . 

Evaluation of the drop-acceleration term a ^; The remaining 

term to be evaluated in equation (4) is the acceleration of the drop 
at the airfoil surface a d . To determine the variation of this term 
with chordwise position, values of a d were calculated from the 
trajectory data by equation (4) for each of the airfoil cases presented 
in tables I through V. The procedure used in making the calculations 
was to compute the value of R/Ry by utilizing values of the 
orthogonal drop— velocity components from tables I through y for 
corresponding values of + and Ry. The term was calculable through 
knowledge of R/Ry and Ry. The terms R/Ry, C^/24, Ry, and V were 
then substituted into equation (4) and solved for a d . The results 
of this procedure are presented in figure 12 for each of the five 
airfoil cases. 

An inspection of figure 12 shows that drop acceleration at the 
surface of the airfoil, like the hodograph of drop velocities for 
tangentially impinging trajectories, can be considered a single rela- 
tion regardless of the combinations of \|r and Ry. How the singular 
nature of the acceleration values arises can be shown by rewriting 
equation (4) as follows: 
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a d 



(5) 


However, since the term (R/Ry)^ is taken to be constant for a 
given position on the surface, equation (5) may be written 


a d = (const) ^C d 


(6) 


where the constant depends on the chordwise position and the value 
of (R/Ry) at that position. Thus, according to equation (6), it 


is apparent that, if the product of if and C d remains constant for 
various values of Ry at a given chordwise position, then the value 
of a d also will remain constant. Hence, to show the constancy of 
a d at various chordwise positions for various Ry values, products 
of if and C d were calculated for every combination of if and Ry 
presented in tables I through V for the five airfoils. The procedure 
used was to select a value of Ry and to derive C d by using the 
value of (R/Ry)^ noted on the hodograph for the s/c position 
corresponding to the particular Ry value chosen. The products of 
\|iC d were compared with other >l<C d products corresponding to the 
same physical point on the airfoil; but these latter products were 
calculated by using values of if obtained from a plot of the tabu- 
lated values of t and s/c for various values of Ry. This 


procedure was adopted in order to obtain products corresponding to 
farthest positions of impingement not given directly in tables I 
through V. Comparisons of the products are made in tables VII 
through XI. These comparisons show that, for a given s/c position, 
the product of \Jr and C d generally is of the same order of magni- 
tude for a wide range of \|r and Ry values, particularly at the 
lower values of s/c. Therefore, the assumption that a d is con- 
stant for a particular chordwise position seems fairly well justified. 


The problem of evaluating drop acceleration now has been 
reduced to finding, for a number of s/c positions on the airfoil, 
a value of a d which corresponds to each s/c value chosen. 

In appr oxima ting the drop acceleration at a point where the 
drop trajectory is tangent to the airfoil surface, several pro- 
cedures are possible, as was the case with the term R/Ry. Of the 
various procedures investigated, the one which will be presented 
herein is considered most acceptable because the resultant accuracy 
is co mme nsurate with that produced by the most accurate procedure 
presented for obtaining R/Ry. In addition, the procedure is simple 
in application. For this procedure, the approximation is made that 
the tangential acceleration of a drop at a given point on the surface 
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is the same as the acceleration of the air along the airfoil surface 
at the same point. The equation used to express the drop acceler- 
ation in terms of air. velocity at the airfoil surface is: 


a d = 



a(u a /v) 

d( s/c) 


(7) 


The velocity gradient term in equation (7) can be evaluated simply 
by plotting U a /v against s/c, and obtaining the slope of the 
resulting curve at the desired s/c positions. 


Results obtained by using equation (7) to approximate values 
of aa for the five airfoil cases are shown in figure 12 where 
the calculated point values are denoted by square symbols. It will 
be noted that equation (7) appears to provide good results for the 
symmetrical Joukovski airfoil case (fig. 12(a)), and over most of 
the lower surface of each of the other four airfoil cases (figs. 12(b), 
12(c), 12(d), and 12(e)). For these latter four airfoil cases, 
equation (7) provides upper— surface drop-acceleration values toward 
the leading edge of the airfoil which appear in good agreement with 
the trajectory data; but farther aft, the ability of equation (7) 
to predict appropriate values diminishes appreciably, particularly 
for the Joukowski and NACA 65 2 -015 airfoils at 4° angle of attack 
(figs. 12(c) and 12(e)). In these instances, the inability of 
equation (7) to yield drop acceleration values fairly feu: aft on 
the airfoil surface apparently is because the drops impinging in 
this region have sufficiently large inertia so as not to respond to 
the very rapid changes in surface— air velocities prevailing near 
the position of maximum air velocity. Except quite near the leading 
edge, the trajectories are fairly straight, indicating that the 
impinging drops do not respond appreciably to the vertical components 
of air velocity. Thus, it seemed that a suitable representation of 
drop acceleration for the case of an airfoil with a velocity peak 
located near the leading edge might be obtained by assuming that the 


Only the tangential component of drop acceleration needs to be 
approximated since the normal component of drop acceleration can 
be shown to be equal to zero at the point of drop impingement 
for tangential trajectories. That the normal acceleration of the 
drop is zero at the point of impingement can be shown by writing 
the equations expressing dynamic equilibrium of a drop. The terms 
involving the drop and air velocities are resolved normally and 
tangentially. A substitution of the boundary conditions at the 
point of impingement shows that the normal acceleration must equal 
zero at this point. 
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drop acceleration was influenced mostly by the x components of air 
velocity. In equation (7)> then U a /V would be replaced by u a /V 
so that 


a 


_ u^ dCug/Y) 

d " V d( s/c) 


( 8 ) 


Equation (8) was used to calculate additional drop-acceleration 
values for the upper surfaces of the airfoils presented in figure 12 
and the results obtained are compared in these figures with drop 
acceleration values given by the trajectory data. The points calcu- 
lated by equation (8) are shown as triangles in figure 12, and 
represent a good approximation of the upper surface a^ values 
fairly far aft on these surfaces. 


The question arises as to whether it would be possible in the 
general case , When the differential analyzer data points shown in 
figure 12 were not present, to detect the inadequacy of equation (7) 
to represent the correct values of a^. In this connection it should 
be noted that the point on the curve for a^ = 0 can always be 
selected because this point corresponds to the chordwise position of 
the tangentially impinging straight-line trajectory having the 
maximum s/c intercept. This particular trajectory can always be 
established by constructing a line tangent to the upper surface of 
the airfoil and parallel to the free— stream velocity. With this 
point located it would be quite evident that the square data po in ts 
(equation (7)) in figure 12(e), for instance, could not represent 
the correct curve, and that equation (8) should be used to obtain 
additional points to provide -a reasonable guide. The technique 
of fairing a curve, through calculated drop-acceleration values, 
to a value of zero acceleration at the extreme position of drop 
impingement can be used also for the lower surface of airfoils; but 
the technique can be applied only in those instances where the 
airfoil angle of attack is sufficiently low to permit a straight- 
line trajectory to impinge tangentially at some point on the airfoil 
surface . 

S ummar y of calculation of scale modulus tfr : The two preceding 

subsections have shown how the terms in equation (4) can be evaluated 
to obtain values of \|r for selected Ry values at chosen positions 
on the airfoil surface. However, a special procedure for evaluating 
if at the stagnation point is necessary, since equation (4) cannot 
be used to evaluate the scale modulus at or very near the stagnation 
point. This procedure is more suitably discussed in connection with 
the section on rate of impingement which follows. 
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For a further verification of the calculation of scale modulus, 
appendix A includes a comparison of results of empirical calculations 
of farthest position of impingement with results obtained from formal 
solutions of the trajectory equations for five airfoils. 

Generalization in regard to rate of impingement .— Another 
quantity of interest to the designer of an aircraft thermal— ice- 
prevention system is rate of drop impingement on an airfoil. An 
expression for weight rate of drop impingement per unit length of 
span, according to reference 7 > is given by 


M s = 3600Vm£y o » 
t 


(9) 


A7o+ 


In order to evaluate the rate of impingement Mg, the term 
must be known. When methods like those of references 3> and 
6 are employed, * can be determined directly from the calculated 

trajectories which impinge tangentially upon the airfoil. For an 
empirical procedure where no trajectories are available, however, 
determination of A 7 0 ' mus ' t be based upon quantities which are 

known. Preceding sections have shown that (s/c)„ and ( 3 / 0)7 can 

L t 

be established as a function of + for various values of Eyj 
hence, the airfoil ordinates corresponding to the farthest position 
of drop impingement on the upper and lower surfaces y^ and y, 

also can be ascertained as a function of ijr for various values of 
By. Because values of y^ and y , can be obtained readily for a 

wide range of \|r and Ry values, the data were examined for a 
relationship involving Ay 0 ^_ (for small angles of attack, £y 0 

is approximately equal to Ay 0 •) and the quantity y„ — y^ which 

t u t t 

will be called Ay+* In this regard, && was compared to £y, 

U O-f; t 

over the range of \|r and By values presented in tables I 
V for the five airfoil cases. 


through 

The results are shown in figure 13. 


An inspection of figure 13 shows that, for each airfoil, the 
ratio of bo Ay t can be considered linear with respect to the 

log of the scale modulus + for various Ry values. The linearity 
exists for values of < 0.8 for the Joukowski airfoils, 

and a value of (£y 0 /Ay).f. < 0.9 for the NACA 6^-015 airfoil; but 
this linearity appears to be characteristic only of airfoils as 
cylinder data from reference 6 , when plotted in the same manner, do 
not show this property. In figures 13(c) and (e), the ratio notice- 
ably exceeds unity, the reason being that the projected frontal 
height h (equal to ^o ^. 1 I° r bhe limiting case of straight-line 

trajectories) exceeds the airfoil maximum thickness when the airfoil 
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angle of attack is other than zero. Of special interest in figure 13, 
however, is the fact that the ratio (Ay 0 /Ay)^. must "become zero at 
some particular value of + for a given value of Ry. This 
"critical" value of \|r can he calculated from an aerodynamic prop- 
erty of the airfoil. According to reference 6, for symmetrical bodies 
at zero angle of attack, the critical value of i|r (i.e., values 
greater than that for which no drops impinge on the body) is given by 


* = 4 Ett ^(ua/v) 

v cr 


$=0 


( 10 ) 


For symmetrical bodies at an attitude other than 0°, or for unsyrn- 
metrical bodies at an arbitrary attitude, the same form of equation 
(.10) applies, but with the notation slightly altered, thus 


t 


cr 


toy 


d(Pa/y) 

dS 


$ =0 


(H) 


This change is made because the small drop which just impinges at 
the stagnation point of the airfoil closely follows the stagnation 
streamline which, in the general case, is not a straight line. For 
simplicity, though, equation (ll) shall be written 

* cr = (12) 

In order to use equation (12), the problem presents itself of 
assig nin g a value of G for the case of an aribtrary airfoil. 
Fortunately, the problem is relieved by the fact that the answers 
sought, s/c and E, are effected only in a minor way by variations 
in G. 4 Thus, it was believed that for determining G the airfoil 
could be replaced by a shape more amenable to calculation. It was 
assumed that a symmetrical Joukowski airfoil would be representative 
of that type section having maximum thickness fairly well forward, 
and an ellipse representative of that type section having maximum 
thickness well aft (low-drag airfoils). Since the major factors 
influencing the value of G are thickness and angle of attack, 

4 Calculations have shown that negligible changes in s/c and E 
incurred for a change in G as large as 10 percent. 
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calculations of G were made for a symmetrical Joukowski airfoil 
and an ellipse of different thickness and lift coefficients. 5 The 
results of these calculations are presented in figure l4. The data 
in 14(a) are intended for use with airfoils resembling Joukowski air- 
foils, and may he used directly. The data in figure 14(b) are 
intended for use with low-drag profiles; however, it is first neces- 
sary to establish an "equivalent ellipse’' thickness ratio for the 
low-drag section being used. An equivalent ellipse is defined for 
the purposes of figure 14(b) as an ellipse having its leading-edge 
radius equal to the leading-edge radius of the airfoil, and a thick- 
ness equal to the airfoil m ax Inr um thickness. The major axis of the 
e -4IdP se is thus established, and the ellipse thickness ratio can be 
computed. An equation expressing the thickness ratio of the equiv- 
alent ellipse in terms of the airfoil leading— edge radius and thick- 
ness ratio is: 



(13) 

_ , W ^ t ?, the aid of figure 1] ^ the value of * for a large number 
of airfoils can be estimated for any Ky value in accordance with 
equation (12). Not only does this value correspond to the condition 
of zero rate of impingement, but it corresponds to the c ond ition 
of zero area of impingement. Hence, the critical value of f can be 
used for obtaining an additional point for area of impingement 
computations, and this value will correspond to the s/c value at 
the stagnation point. 


While the condition of no drops impinging on the airfoil surface 
yields one point on the curves, (Zy Q /Ay) t versus log + # at 

least one more point is required for each value of R v in' order to 
establish the linear relationships observed in figure 13. To locate 
a second point on an isopleth of By, it is desirable to determine 
a value of * corresponding to a chosen value of (Ay,//y) + somewhat 
less than unity. The reason for this specification is to procure a 

rjf 53 ° f , <V^>t used to establish the linear 

^ i0nShiPS A et T een (Ay o/^) t aa* log *, for isopleths of R. 
which were noted in figure 13. V 


5 No account is taken of the effect of a cambered profile on the 
velocity gradient G. The reason for this is that tests with an 
electrolytic tank have shown that the effects of camber are very 
small in comparison with the effects of thickness, and calculations 
have shown that only large variations in G are important in 
affecting the values of s/c and E. 
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In developing a procedure for determining what value of If 
is associated with a specified value of (Ay 0 /Ay) t on an isopleth of 
Ry, the data from the five airfoil cases were examined for values 
of some parameter, related to (s/c)^ and (s/c)j^_, which could 

be used to fix the value; of ♦ . The parameter used to supply the 
necessary values was the efficiency of drop impingement E. The rela- 
tionship between E and (Ay 0 /Ay) t is given by 



Equation (l4) can be derived by starting from the definition of E 
in terms of the initial drop trajectory ordinates 



At the small angles of attack associated with most flight conditions, 
h in equation (15) can be replaced by and Ayo’ *>y AFo t 

so that 


^o. =Et max (16) 

"C 

Then, if both sides of equation (l6) are divided by Ay^ and the terms 
rearranged, equation (l4) is obtained. 

The trajectory data for the five airfoil cases provided, for dif- 
ferent values of Rv, relatively constant values of E correspond- 
ing to a value of 6 (Ay o My)^. = 0.8. These efficiency values were 

6 The procedure utilized was to determine from curves of (Ay 0 /Ay)t, ™ 
as a function of If (fig. 13), the value of at which 
(Ay 0 /Ay).f. =0.8 for different values of Ry. Then, data from tables I 
through V were used to establish curves of E , as a function of t 
for the same values of Ry. On the efficiency curves, the value of E 
corresponding to - (Ay 0 /Ay)^. =0.8 for a particular value of Ry 
could be determined by locating, for the same Ry value, the value of 
If ' which was established from curves of figure 13 to correspond to 

(Ay 0 /Ay)t = °* 8 * 




22 


NACA TN 2476 


used to obtain an average efficiency value for each airfoil case. 
Then, by using equation (l4) , an average value of Ay t /t Tnflx could 
be computed for each airfoil case by using the average efficiency 
values and a value of (A y 0 /Ay) t = 0-8- The results are presented 
in the following table: 


Case 

number 

Efficiency of 

impingement 

, E 

(percent) 


Ay t 

t 

max 


% 

Average 
value for 
each case 

16 

32 

64 

128 

256 

512 

1024 

1 

2048 

I 

— 

77.5 

— 

77.0 


75.5 

— 

78.0 

77.0 

0.963 

II 

76.0 

— 

74.5 

' 

75.5 

— 

77.5 

— 

75.8 

.948 

III 

72.0 

— 

73.0 

— 

72.0 

— 

70.5 

— 

71.8 

.898 

IV 

77.0 

— 

82.0 

— 

86.0 

— 

82.0 

— 

81.7 

1.022 

V 

55.0 

— 

59.0 

— 

56-. 5 

— 

55.0 

— 

56.4 

.705 


The values of Ayt/’ t max "tabulated in the preceding table exhibit 
some variation between airfoil cases, and figure 15 is presented to 
show this variation when AFtAmax ia assumed to be a function only of 

angle of attack. In figure 15, a curve is presented for the Joukowski 
airfoils, but not for the NACA 652~OT5 airfoil because the trajectory 
data provides only one point. Until additional information is avail- 
able on the' variation of Aft Amax of attack, it is recom- 
mended that the variation shown in figure 15 be used. It is possible, 

then, to determine, for a given value of Ry, an approximate value of* 

\|r at which (/2y 0 /Ay)^ = 0.8. The procedure which may be used for 

determining this value of \| r is shown graphically in figure 16. From 
curves of (s/c)^ and (s/c) lt as a function of log i for a 

specified value of Ry (fig. 16(a)), curves of y and y^ as a 

function of log \|r are established for the same value of Ry 

(fig. 16(b)). For the relation shown by figure 16(b), there is a 

value of Ayt AmflY which ia the same as the one chosen from the rela- 

tion in figure 15 corresponding to the airfoil angle of attack. This 
particular value of Ay-tAmax corresponds to the ^ value at which 

(Ay o /Ay) t = 0.8 for the particular R^ value chosen (fig. 16 (c)), 
and the second point on an isopleth of R y for (Ay 0 /Ay)t as a 
function of log i|r is thereby determined. 
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The previous discussion has shown how values of Ayo-^ ®&y be 
obtained for various and Ry values. However, in the design 
of a thermal ice— protection system, by the method discussed in 
reference 1, it is sometimes more convenient to determine the rate of 
water-drop impingement by using the airfoil collection efficiency E 
rather than by using the term Ay 0t ' . In such circumstances, 
equation (9) becomes 

M s = 3600 Vm E t^ 

wherein E is given by equation (l4). When equation (l4) is used and 
the angle of attack is other than zero, the limit efficiency value 
corresponding to straight-line trajectories will be greater than unity 
because h usually is somewhat greater than t ^y . 

Appendix A includes a comparison of results of calculations of 
efficiency of impingement by the empirical method of this report with 
results obtained from formal solutions of the trajectory equations for 
the five airfoil cases. 

Generalization in regard to distribution of impingement .— Of 
secondary importance in the design of heated wings is the distribution 
of water-drop impingement over the length of interception along the 
airfoil surface. However, knowledge of distribution of water drops 
over an airfoil is sometimes quite useful and, therefore, a discussion 
of the subject is believed worthwhile. 

An examination of the trajectory data did not reveal any direct 
way to obtain an empirical functional relation between impingement 
distribution, scale modulus, and free— stream-drop Reynolds number. 
Therefore, a graphical construction was resorted to in order to approx- 
imate the distribution of drop impingement over an airfoil surface. 

That a graphical procedure could be used in this manner was found by 
examining plots of the concentration factor 7 C as a function of s/c 
for various combinations of t and Ry. Typical variations of these 
factors are presented in figure 17 calculated from data for the five 
airfoil cases. For each of the airfoil cases presented in figure 17, 
two curves represented by solid lines are shown. One curve is typical 
for combinations of ^ and Ry corresponding to curved trajectories, 
and the other curve is typical for the combination of \|r and Ry 
corresponding to straight-line trajectories (t=0, value of Ry 
arbitrary) . The curve for i|r=0 is obtained by drawing a number of 
straight-line trajectories to the airfoil to obtain values of the 
concentration factor, C = dA 0 /dA s , and represents the locus of 
msTiimim possible values of C. This curve, which will be referred to 
as a limit curve, can always be obtained for a given airfoil because 

7 The use of the concentration factor C in the computation of heat 
requirement due to drop impingement is discussed in reference 1. 
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straight— line trajectories can. always be reproduced, but the curves for 
values of C less than maximum cannot be obtained because the shapes 
of the curved trajectories corresponding to these values c ann ot be 
determined. Because of the shape of the C distribution curves noted 
for the five airfoil cases in figure 17, a triangular distribution 
has been found, useful to represent approximately the distribution 
curves for curved trajectories. If a triangular distribution is 
assumed, as shown by a dashed line for each of the airfoil cases in 
figure 17; the following equation may be written for the ma-sriTmrm value 
of C on the triangular distribution: 



1 

2 C max 



+ 



from which 



(IT) 


The denominator of equation (17) can be evaluated from the curves of 
farthest positions of tangential drop impingement. The numerator may 
be evaluated from the equation 


so that 



( 18 ) 


-'max 



(19) 


That the relation expressed in equation (18) exists can be shown 
by recalling from equations (9) and (15) that the total weight rate of 
drop impingement on an airfoil may be written as 


M s - = 3600 Vm Eh (20) 

Equation (20) may also be written as the summation, between the farthest 
positions of impingement, of the local weight rates of impingement per 
unit surface area, thus . . 
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T Ut Mds- (2D 

Je h 

in equation (21) can be expressed by 
Ma = 3600 VmC (22) ' 

where 

C = dA 0 /dA s (23) 

By substitution of the value of Mq from equation (22) into 
equation (21), the result is obtained that 


M 0 = 


According to reference 1, Ma 


M g = 3600 Vm 


r* 


Cds 


(2k) 


Then, by comparing equations (20) and ( 2k ), equation (l8) can be 
obtained. 

Use of equation (19) permits the calculation of for a 

triangular distribution because all the terms of the right-hand member 
are known. The value of is taken to be on a line connecting the 

points C=1.0, s/c=0, and C=0, and s/c for the stagnation point. 

The reason for this assumption is that the very large drops are very 
little deviated from their paths, and the very small drops tend to 
follow the stagnation streamline. 

The value of C mnY obtained from equation ( 19 ) always will be low. 
However, if the triangular approximation is altered to correspond more 
nearly to the shape of the limit curve for the C values, while 
retaining the enclosed area the same as the triangular area, more 
accurate concentration-factor values can be obtained. The altering 
of the curve is an attempt to establish the locus of concentration^ 
factor values as it would be given by data obtained for calculated 
trajectories. , . - 

Generalization in regard to calculating area and efficiency of 
impingement at high subsonic speeds .— With the advent of airplanes in 
the 500 to 600 miles per hour speed range, considerable interest has 
been shown in the effects which compressibility may have on area and 
rate of drop impingement. 
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The differential analyzer computations for the five airfoil cases 
of this study were made to encompass airspeeds up to about 350 miles 
per hour on the basis that any effects of compressibility could be 
neglected. Consequently, no corrections to the velocity field were 
introduced for any of the airfoils to account for compressibility. 
However, for higher speeds, the effects of compressibility would be 
expected to become more pronounced and some method of taking compressi- 
bility into account should be considered. 


The effects of compressibility : In compressible flow below 

the critical Mach number, streamlines at a large distance forward of 
the airfoil have nearly the same shape as for incompressible flow. 
Wearer the airfoil, deviations of the compressible flow streamlines 
from the incompressible flow streamlines are greater (reference 8). 
Thus, in compressible flow, drops may be considered to start at 
infinity at the same initial conditions as they would for incompres- 
sible flow, but terminate on the airfoil surface at velocities 
governed by the local compressible— flow surface-air velocities. 

These velocities are, in general, higher- than corresponding incompres- 
sible-flow velocities; therefore, compressibility may have some 
effect on farthest position and efficiency of drop impingement. 


If it is presumed that drops in compressible flow react to a 
change in surface-air velocities (due to changes in flow field 
arising from compressibility) in much the same manner as drops react 
to changes in flow field caused by a change in airfoil basic thick- 
ness, load distribution, or angle of attack, the same general 
procedure used to calculate farthest position and rate of imp ing ement 
for incompressible flow can be applied to the calculations for 
compressible flow. 

To obtain the air hodograph for compressible flow when experi- 
mental data are not available, the equation 



L 


1 + 


l-(l + 0.7025M 2 Cp M ) 
0.2025M 2 


O .2883 


1/2 


(25) 


(taken from reference 9) “ay be used to evaluate the magnitude of the 
velocity vectors. The pressure coefficient, Cp M , in equation (25) 
can be calculated using the Prandtl-Glauert equation 
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where . 

P = A /l-M 2 

and Gp is the pressure coefficient at a point on the airfoil surface 
when the flow is incompressible. The pressure coefficient Cp is 
related to the surface— air velocity at the same point on the surface 
by the equation 

/u a \ 2 

C P “ 1_ (“ ) (26) . 

Values of U a /v used in equation (26) to calculate various values of 
Cp are the same as those used to establish the air holograph for 
incompressible flow. 

Effect of compressibility on the farthest position and effi- 
ciency of drop impingement : Because of incompressible— flow holo- 

graph is altered to another shape for compressible flow, the drop 
hodograph also will be modified, utilizing the procedure developed 
for the incompressible case. The changes in the drop and air holo- 
graphs in turn alter the values of a d and (R/Ry) at points along 

the airfoil surface. According to equation (4), if changes in a d 

and the product (B/Ry)^ (C^b/ 24) are not proportionate, the scale 
modulus will be altered for given ( s/c) t and Ry values. To show 
the effect of compressibility on the scale modulus, some calculations 
of farthest position and efficiency of drop impingement were made for 
the 15-percent— thick Joukowski airfoil at zero lift and M=0.6 (which 
is near the critical value) . A value of Ry=225 was chosen for the 
calculations (430 mph, 15,000 ft, 0° F, and drop diameter of 25 
microns). The results of these computations are compared in 
figure 18 with empirical results for incompressible flow. Figure 
18 (a) , which presents the farthest position of drop impingement as 
a function of scale modulus, shows that the scale modulus has been 
increased slightly over most of the range in ( s/c) t for the isopleth 
of By chosen. This result is equivalent to the farthest position 
of impingement being increased from the value of (s/c)t corre- 
sponding to the same i|r . and Ry value in incompressible flow. For 
the conditions chosen, the increase in farthest position of. impinge- 
ment is small and, for most practical purposes, may be neglected. 
Figure 18(b) shows that compressibility also slightly increases 
impingement efficiency. This result can be explained by equation (l4) 
which indicates that increases in E are effected by increases in 
either the value of (Ay 0 /Ay)t or AFt> or b °th. Increases in the 
values of both of these terms, for given . ifr and By values, occur 
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for increased values of farthest position of impingement, and hence 
impingement efficiency is increased. The increase in efficiency due 
to compressibility for the same airfoil also is small and; like the 
increases in farthest position of drop impingement, may be neglected 
for most practical purposes. It should be noted, though, that the 
computations for farthest position and efficiency of impingement 
were made for a relatively low Mach number (0.6), and that somewhat 
larger changes in both parameters for given \|r and Ry values may 
be possible for airfoils having appreciably higher critical Mach 
numbers. However, it is suggested that the empirical method as 
applied to compressible flow not be used for airfoil sections having 
high critical Mach numbers (0.9) in view of the fact that no data 
exist for giving validity to the extrapolation of the method of this 
report to compressible flow. 


PROCEDURE FOR USING THE METHOD TO CALCULATE AREA, 
RATE, AND DISTRIBUTION OF WATER-DROP 
IMPINGEMENT ON AN ARBITRARY AIRFOIL 


Previous sections have shown how the empirical method of this 
report may be applied to determine area, rate, and distribution of 
impingement for an arbitrary airfoil in either incompressible or 
compressible flow. The procedure to be followed in using the method 
will now be explained by means of an example of the computations in 
the case of an NACA 23015 airfoil, in incompressible flow, at c ? =0.5. 
The detailed computation procedure in this example is arranged in a 
form believed to be most readily usable; and the various phases of the . 
computations are presented in a step-wise manner, each phase consisting 
of several operations for ease in following the procedure. The 
procedure for compressible flow is nearly the same as. for incompressible 
flow, hence, only the procedure for incompressible flow is illustrated* 
However, where differences in procedure would occur for the case of 
compressible flowj these are noted. 


Area of Impingement 


The procedure for calculating area of impingement consists pri- 
marily in determining values of (s/c)^ and (s/c)^. The following 

steps explain how the empirical relations derived from the trajectory 
data are used to determine these values, and figure 19 incorporates 
necessary accompanying graphical relationships. 
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Step 1: The following curves are needed for use during the 

computation procedure: 

(a) A large-scale plot of the airfoil (fig. 19(a)) 

(b) A plot of s/c versus x for both upper and lower 
surfaces (fig. 19 (b)) 

(c) A plot of k for various x positions (fig. 19(c)) 

(d) Chordwise distribution of incompressible— flow air 
velocities over the airfoil surface (fig. 19 (d)) 

Step 2: Construction of the air hodograph: • 

(a) On a coordinate system using vertical air— velocity 
components as ordinates and horizontal velocity components 
as abscissas, measure values of air velocity as vectors 
emanating from the origin. (These values are obtained from 
the velocity distribution (fig. 19 (d)) established in 

step 1(d)). The slope of the velocity vectors with respect 
to the horizontal axis is determined by the shape— factor 
curve established in step 1(c); Figure 19(e) illustrates 
this step for the NACA 23015 airfoil. 

(b) Connect the end points of all vectors noting that the 
envelope curve (the hodograph) passes through the origin 
for the physical coordinates corresponding to the stagnation 
point. Refer to figure 19(e) for illustration. 

(c) For compressible flow over a straight wing below the 
critical Mach number > compute the velocities for the air 
hodograph from 

1— ( 1+0 . 7025M 2 Cp ) 0,2883 

_ r M 

0.2025M 2 


Step 3: Construction of drop hodograph: 

(a) From the hodograph established in step 2(b) 
(incompressible flow) or step 2(c) (compressible flow), 
note the maximum value of Vertical air— velocity component 
Va max/ V f°r. ..both upper and lower surfaces. For the NACA 

23015 airfoil, the air hodograph shows that v Q /v for 

a max / 
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the upper .surface is 0.915 and. for the lower surface , O.39. 

(b) Using figure 10, determine the value of (va^^/y) — y /V) 

for both surfaces. A subtraction of fcajnax^ “ ( v d x / v ) from 

v a m AY /V will yield the vertical location of v d ffla ^?^* Project 

v^ax/V horizontally for each surface until it intersects the ray 

connecting the origin and v a /V. In this manner, the value 

max 

of (R/R v ) t corresponding to Vg ^y /V is determined graphically. 

In the case of the NACA 23015 airfoil, (v„ /V) -(v^ /v) for 

v d max' ' v a max' ' 

the upper surface is 0.245 and for the. lower surface, 0.068. 

These values are measured from the air hodograph (fig. 19(f)) 
to provide values of (R/Ry) corresponding to upper and lower 
surface values of v„ /v. 


Tnax 


(c) Establish the vector passing through the origin and inclined 
with the horizontal axis at an angle equal to the airfoil angle 
of attack. The length of this vector is made equal to unity 
(fig. 19(f)). . • 


(d) For both airfoil surfaces, establish an approximate 
tangential— trajectory hodograph which follows the equation 



WL 

(u a A) v , 


(27) 


The hodograph determined by equation (27) for the upper 
surface of the NACA 23OI5 airfoil is designated as curve A 
in figure 19(g) and as curve B for the lower surface. When 
equation (27) provides, for either the upper or lower surface, 
a drop hodograph' which fails to pass through the point 

= cos <x, v d /V = sin a (point X in fig. 19(g)), a curve 
should be-faired from point X to connect tangentially with the 
approximate • drop hodograph as determined by equation (27). 

In the example case, neither the upper nor lower surface 
drop hodographs passes through point X, and therefore curves 
are faired from this point to the approximate drop hodographs. 
The faired curves are labeled C and D in figure 19(g) for the 
upper and lower surfaces, respectively. For the upper surface, 
the approximate drop hodograph almost ; passes through point X, 
so little difficulty is experienced in fairing a connecting 
ourve. However, for the lower surface, the' approximate drop 
hodograph deviates considerably from point X leaving a larger 
latitude for fairing a curve from this point than was the 
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case for the upper surface. In this instance, the contour 
of the air hodograph was followed as a guide in fairing a 
curve from point X to the approximate drop hodograph. The 
combination of curves A and C for the upper surface, and 
curves B and D for the lower surface, comprise the resultant 
tangential— trajectory hodograph for each surface. 


(e) Convert the x values used in the air hodograph of 
figure 19(e) to s/c values by utilizing the; relationship 
in figure 19(b) between x and s/c. This step is shown 
also in figure 19(g). 


(f) Measure the values of (s/Ey) 


t 


as determined by the 


drop and air hodographs in figure 19(g) and plot these 
values against their corresponding s/c positions. This 
step is performed in figure 19(h) for the upper and lower 
surfaces of the NACA 23OI5 airfoil. 


Step 4 : Calculation of drop-acceleration values at the air- 

foil surface. 

(a) Use the air— velocity distribution over the airfoil 
surface (fig. 19 (d)) and the relation between x and s/c 
(fig. 19 (b)) to establish the relation required between 
U a /V and s/c for both upper and lower surfaces. 

Equation (7) then can be used to calculate values for a^. 
Values obtained for a^ by equation (7) are plotted 
against their respective surface positions s/c, and this 
is done in figure 19(i) for the upper and lower surfaces 
of the NACA 23015 airfoil. 

(b) Include- in the plot of a^ versus s/c (fig. 19(i)) 
the value of s/c at which a^ is zero on each surface. 

The value of s/c at which a^ is zero corresponds to 
the extreme position of tangential drop impingement on 
the airfoil surface. For each surface, the extreme 
position of drop impingement can be obtained by drawing 

two straight-line trajectories in the free-etream direction, 
one tangent to the airfoil upper surface and the other 
tangent to the lower surface. The point of tangency for 
each straight-line trajectory defines the s/c . value 
corresponding to an a^ value of zero. For the NACA 23015 
airfoil at a lift coefficient of 0.5j the extreme position 
of impingement on the upper surface is at s/c=23.0 and 
on the lower surface, at s/c=53.0. 
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(c) For each surface, fair a curve through the calculated 
values of a^ to the value of s/c found in step 4 (b) 
corresponding to a d of zero. This step is also performed 
in figure 19(i). On the upper surface, values of a^ given 
by equation ( 7 ) were considered in error over a small range 
of s/c near the peak in the velocity distribution curve 
(fig. 19 (d)), and consequently little weight was given the 
points in this region when extending the faired curve to 
a^=0 at the extreme position of drop impingement'. A few 
additional points can be calculated by equation ( 8 ) to aid 
in the fairing of the drop-acceleration curve, but since in 
this example the extreme position of impingement is fairly 
far forward (b/c = 23 .O) and equation ( 7 ) will yield posi- 
tive a^ values for s/c values as far aft as 11 percent, 
use of equation ( 8 ) was deemed unnecessary. 

Step 5: Computation of scale moduli corresponding to arbi- 

trarily selected values of farthest position of impingement. 

(a) Using drop-acceleration values obtained from figure 
19 (i) of step 4 and (R/Ry)^ values from figure 19(h) of 

step 3, calculate value,s of scale modulus from equation 
(4). The computations are made for arbitrarily selected 
values of Ry and s/c. Values of C d R/24 required for 
solving equation (4) can be obtained from interpolation of 
the data presented in table VI, or preferably from a semi- 
logarithmic plot of the data. A small range of tabulated 
C^R/24 values, including the values involved in the compu- 
tation of V for the NACA 23 OI 5 airfoil is shown graphically 
in figure 19(j). The computations of \|r for the NACA 23 OI 5 
airfoil are shown in the following table: 



NACA TN 2476 


33 


CALCULATION OF SCALE MODULI FOR 
CHOSEN AREAS , OF IMPINGEMENT AND DIFFERENT R y VALUES 

Ry=10 

s/c 

a d 

a d% 

(R/%) t 

R 

C^R 

24 



UPPER SURFACE 

5 

10 

15 

20 

8.6 

2.8 

1.0 

.2 

86 

28 

10 

2 



0.365 

.405 

.412 

.405 

3.65 

4.05 

4.12 

4.05 

1.42 

1.45 

1.46 
1.45 

.518 
.586 
.600 
' .586 

166 

48 

17 

4 

LOWER SURFACE 

5 

10 

20 

30 

ho 

45 

2.7 

1.0 

.6 

.35 

.17 

.10 

27 

10 

6 

3.5 

1.7 

1.0 

0.10 

.12 

.10 

.085 

.072 

.070 

1.0 

1.2 

1.0 

.85 

.72 

.70 

1.17 

1.20 

1.17 

1.14 

1.13 

1.12 

0.117 
.144 
.117 
. .097 
.081 
.078 

231 

70 

51 

36 

21 

13 

Ry=100 

UPPER SURFACE 

5 

10 

15 

20 

8.6 

2.8 

1.0 

.2 

860 

280 

100 

20 

0.365 

.405 

.412 

.405 

36.5 

40.5 
41.2 
40.5 

2.88 

3.00 

3.02 

3.00 

1.05 

1.21 

1.24 

1.21 

819 

230 

80 

16 

LOWER SURFACE 

5 

10 

20 

30 

40 

45 

2.7 
1.0 
.6 
• 35 
.17 
.10 

270 

100 

60 

35 

17 

10 

0.10 

.12 

.10 

.085 

.072 

.070 

10 

12 

10 

8.5 

7.2 . 

7.0 

1.80 

'1.90 

1.80 

1.72 

1.64 

1.63 

0.180 

.228 

.180 

.146 

.118 

.114 

1500 

439 

333 

240 

144 

88 


Step 6: Establishing the curves of farthest position of 

impingement . 

(a) Plot the values of \| 1 calculated in step 5 against 
the corresponding values of s/c as ordinates for isopleths 
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of Ry. Use a semi logarithmic- scale, plotting the ordinates 
arithmetically. This is done for the NACA 23OI5 airfoil in 
figure 19(h). 

(b) On the plots of s/c begun in step 6(a), denote by 
lines drawn at constant s/c values the value of extreme 
position of drop impingement on eaoh upper and lower sur- 
face . For the NACA 23OI5 airfoil, at a lift coefficient 

of 0.5, the extreme positions of impingement have been found 
in step 4(b) as being at s/c=23.0 on the upper surface, 
at s/c=53»0 on the lower surface. 

(c) Calculate from equation (12) the critical— scale modulus 
value for each isopleth of Ry. The value of G required 
to evaluate equation (12) for either incompressible or 
compressible flow can be obtained from figure 14. The value 
of t cr calculated for a given value of Ry corresponds 

to s/c at the stagnation point. Values of ^ cr calculated 
from equation (12) for different values of Ry also should 
be plotted on the graph of farthest position of impingement 
(fig. 19(h)) begun in step 6(a). Since the NACA 23OI5 air- 
foil can be considered to have the same type of profile as 
Jouhowshi airfoils, figure 14(a) is used to obtain a value 
of G for evaluating equation (12). For a value of c^=0.5 
and a thichness ratio of 0.15, figure 14(a) shows the value 
of G to be 37.5. When this value of G is used in 
equation (12), ^cr has the following values: 


Ry 

^cr 

10 

1500 

100 

15000 


(d) For each isopleth of R, fair a curve through the 
values of s/c plotted in steps 6(a) and 6(c) (fig,19 (h)). 
These curves should approach asymptotically, as a limiting 
value, the farthest position of drop impingement for the 
upper and lower airfoil surfaces. Figure 19(h) then may 
be used for design purposes by determining, the values of 
(s/c) u and ( s/c) ^ for particular values of i|r and Ry. 


Rate of Impingement 


The procedure for the determination of the total rate of impinge- 
ment, as has been explained in reference 1, consists of summing the 
rate of water-drop impingement for each of the drop sizes in an assumed 
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drop-size distribution. This is possible for each size of drop by 
use of the equation: 


M s = 3600 EV myinax 


The values of V, m, and y^y are directly obtainable from a knowledge- 
of the nature of the icing conditions and the airfoil shape. The 
procedure for calculating efficiency of impingement consists essentially 
of evaluating equation (l4). The following steps explain how the 
empirical relationships derived from the trajectory data are used to 
solve equation (14), and figure 20 is presented to show the necessary 
accompanying graphical relationships. 

Step 1.— The following relationship should be established 
for use during the computational procedure: s/c as a 

function of y/c for both upper and lower surfaces. This 
relationship is shown for the NACA 23OI5 airfoil in figure 
20(a). 

Step 2.— Establishing (^y Q //iy) t as a function of t for 
isopleths of Ry . 

(a) Using the farthest position of impingement curves 
(fig. 19(k)), establish the airfoil ordinates at points 
of tangential drop impingement as a function of log t 
for isopleths of Ry. This is done in figure 20(b) for 
the NACA 23015 airfoil for Ry values of 10 and 100. 

(b) On the curves established in figure 20(b), establish 
the values of ^ at which (/ty o //ty) t =0 .8. This can be 
done by using figure 15 which provides, for different angles 
of attack, values of Zy t- /t mBY corresponding to the value 
of (Ziy o /Ay)^=0.8. For the NACA 23015 airfoil (assuming 
the NACA 23OI5 airfoil to be similar in profile to the 
Joukowski airfoils), figure 15 gives AytAnmr 0 *^ at 
0=3.6°. This value can be used in figure 20(b) to deter- 
mine the values of \|r corresponding to (£y Q /£tf)^=0 .8* by 
locating the i|r values at which ^yt =0 *9 t max* For tile 
NACA 23015 airfoil, these t values in figure 20(b) are 
taken to be: 



10 

100 


JL 

28 

130 
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(c) Construct a semilogarithmic plot, assigning values 
of scale moduli as abscissas (logarithmic scale) and 
values of (Ay 0 /Ay) t as ordinates (arithmetic scale). 

Then, for each particular chosen value of By, plot the 
values of xjr cr corresponding to a value of (Ay o /Ay)t=0 
and also the values of corresponding to a value of 
(&y o /&y)-t=0.8. This is done in figure 20(c) for the 
NACA 23015 airfoil. In figure 20(c), the values of 'l' cr 
are denoted by circles and the values of t at 
(Z\y o /Ay)^.=0.8 are denoted by squares. 

(d) Using the points plotted in step 2(c), co nne ct with 
a straight line each pair of points having the same . By 
value. This procedure establishes (Ay Q /Ay) t as a function 
of \|r for isopleths of By and is illustrated in figure 20(c) 
for the NACA 23015 airfoil. The curves in figure 20(c) have 
been extrapolated to approach a limit value of (Ay 0 /Ay)^.=1.13 

which corresponds to h/t max =1.09 for straight-line tangential 
trajectories. 

Step 3. -Calculation of impingement efficiency: 

(a) Equation (l4) now is used to calculate impingement effi- 
ciency for various t and By values. To evaluate < 

equation (l4) for a selected value of By, record values of 
(Ay 0 /Ay)t f° r arbitrary values of i|r from the relationship 
(fig. 20(c)) established in. step 2(c). For the same values, 
obtain corresponding Ay t values directly from the curves 
established (fig. 20(b)) in step 2(c). 

(b) Plot the values of E obtained in (a) above against the 
corresponding values of scale modulus in isopleths of By. A 
semilogarithmic scale is recommended, with values of E being 
plotted as ordinates on the arithmetic scale (fig. 20(d)). 

Figure 20(d) then may be used for design purposes by determining 
values of E for particular values of \|r and By. 

Distribution of Impingement 

Distribution of impingement is considered defined, as explained in 
reference 1, when values of the concentration factor C are determined 
over the region of drop impingement. The empirical procedure consists 
first of evaluating distribution of impingement for straight-line 
trajectories. By establishing the distribution of impingement for 
straight— line trajectories, a limit curve (the locus of ttuvsH rmim values 
of C) is obtained which can be used as a guide modifyin g an approximate 
distribution for the particular trajectories of interest. Figure 21 
is presented to illustrate the steps involved for the NACA 23015 airfoil. 
These steps are subsequently described. 
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Step 1.— Impingement distribution on an airfoil caused by 
straight-line trajectories (the limit curve). 

(a) To determine the impingement distribution on an 
airfoil for straight— line trajectories, the procedure is 
to plot the airfoil contour and draw several straight-line 
trajectories, in a streamwise direction, so that they 
impinge upon the airfoil. The two tangential trajectories 
should be included. This procedure is illustrated in 
figure 21(a) for the NACA 23OI5 airfoil. 

(b) For each trajectory appearing in the plot made in (a), 
note the trajectory starting ordinates and the corresponding 
positions of impingement, s/c. Plot these values, using 
position of impingement s/c as abscissas, and fair a 
curve through the points (fig. 21(b)). 

(c) Establish concentration factor values for arbitrarily 
selected values of position of impingement by obtaining 
the slope, at various points, of the curve plotted in (b). 

A plot of these concentration factor values against position 
of impingement values yields the limit curve corresponding 
to straight— line trajectories. Rote that the maximum value 
is always unity at the airfoil leading edge, and the minimum 
values are zero at the extreme position of tangential drop 
impingement . 

Step 2 .—Determining impingement distribution for an arbitrary 
value of scale modulus and free— stream Reynolds number. 

(a) Calculate a. maximum value of concentration factor 
from equation (19). Values of s/c and E needed to 
evaluate this equation for a selected value of Ry and 
various values of i|r may be obtained from efficiency and 
farthest position of impingement curves already established. 
For the RACA 23OI5 airfoil, the value of C maY becomes 


^max 


Eh _ (0.75) (0-164) 
( a av )t l/2(0.135+0A5) 


for a value of =100 and a value of Ry=100 . 

(b) Plot the value of obtained in (a) against 

position of impingement, s/c. Assume that lies 

on a straight line connecting the points 0=1.0, s/c=0, 
and 0=0 for s/c at the stagnation point. See figure 21(c)- 
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(c) Using the value of obtained in step (b), 

form a triangular distribution of drop impingement 
by drawing straight lines from the point 

to the points (a/c)^, C=0 and (s/c)j t , C=0. 

This step is also illustrated in figure 21(c). 

(d) Modify the triangular distribution established in 
(c) to conform to. the general shape of the distri- 
bution curve established in step l(c) for straight- 
line trajectories. In performing this modification, 
the value of C^y is usually increased such that 
the new curve contains the same area as did the 
triangular distribution. The modified curve is shown 
in figure 21(c) for a particular combination of ^ 
and Ry. For other combinations of \|r and Ry that 
may be needed for design purposes, other figures 
similar to figure 21(c) would have to be constructed. 


EVALUATION OF THE EMPIRICAL METHOD 
DEVELOPED IN THIS REPORT 


The degree to which the final values of farthest position and 
efficiency of drop impingement as calculated herein depend upon the 
accuracy of determination of the intermediate quantities (R/Ry) t , 
a,, and G was investigated by determining the effect of arbitrarily 
altering these three quantities a given percentage. With this pro- 
cedure, the effect on farthest position and efficiency of impingement 
can be appraised for the selected changes in the three variables; 
also, . some concept can be gained of the approximate error in the 
empirical method itself. When computations were made for the 
15 — P er cent— ■ thick symmetrical Joukowski airfoil at a = 4°, end the 
values of (R/Ry)±> a d , and G were altered by ±10 percent in all 

possible combinations, it was found for farthest position of impinge- 
ment that in no case was changing G significant. The combination 
of positive and negative changes providing the largest change in f 
resulted in a change in s/c of about 2-percent chord over most of 
the range in values of . The empirical method contributed an 
additional change of only about l/2-percent chord. For efficiency of 
impingement, the effect of a change in the term G alone was to make 
a change in efficiency of about 0.5 percent; the combination of positive 
and negative changes in (R/Ry) and a d providing maximum change 
in made a change in efficiency of about 3 percent over most of the 
range in \|r values. As compared with these changes, the approximations 
of the empirical method led to efficiency of impingement values which 
differed from the differential analyzer .values by about — 15 percent. 
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While the foregoing values will not he necessarily representative 
for all other airfoils, they nevertheless seem to indicate the order 
of magnitude of error in area and efficiency of impingement to be 
expected when the error in the terms (R/Ry)^, a^, and G can be kept 
within ±10 percent. Whether this sort of accuracy always can be 
realized by using the procedures suggested in this report can be 
ascertained only as more water-drop trajectory data become available. 


CONCLUDING REMARKS 


Results of water-drop trajectory data obtained from a differential 
analyzer have provided a basis for the initial development of a 
general method for empirically determining area, rate, and distribution 
of water-drop impingement on airfoil sections of arbitrary profile. The 
method, as applied to the incompressible-flow regime, is more firmly 
established for airfoils resembling the Joukowski airfoils investigated 
than for low-drag airfoils, since the basic data were obtained for four 
Joukowski airfoil cases and only one low-drag section. Because the 
differential analyzer data were obtained only for cases involving 
incompressible flow, the method as applied to compressible flow will not 
provide area and rate-of— impingement data having as much validity as 
the data for incompressible flow. There is no doubt that further water- 
drop trajectory data are needed, particularly for thin airfoils (order 
of 5 — percent thick) at high speeds, and airfoils at high angle of attack 
(in the neighborhood of 12°). Whether these new data would revise the 
procedure presented herein, replace it, or substantiate it remains to 
be seen. Until such data are available, however, the method of this 
report should provide more complete and accurate information on the 
rate, area, and distribution of water-drop impingement on an arbitrary 
airfoil than aAy other known empirical means. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., May 8 , 1951. 
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APPENDIX A 


COMPARISON OF RESULTS OF EMPIRICAL CALCULATIONS OF 
FARTHEST POSITION AND EFFICIENCY OF IMPINGEMENT 
WITH RESULTS OBTAINED FROM FORMAL SOLUTIONS OF 
TRAJECTORY EQUATIONS FOR SIX AIRFOILS 


After a new design method has been developed based on a limited 
amount of data, it is desirable, but not always possible, to evaluate 
the accuracy of the method by applying it to cases where checks on the 
method exist. Unfortunately, in this investigation no experimental data 
exist and. very little calculated data are available on water-drop trajec- 
tories beyond that reported herein. However, there is some justification 
for checking the results of this method against the more accurate calcu- 
lations for the five airfoil cases since these data (except for distri- 
bution of impingement) were used only to establish a procedure which, 
when applied, would not necessarily yield the same values for area and 
percent impingement as the data. Consequently, the empirical relation- 
ships just .established were applied to the five airfoil cases, and 
values of farthest position and efficiency of impingement were then 
obtained and compared with the results of the differential analyzer 
study. These comparisons for farthest position and efficiency of im- 
pingement are shown in figures 22 and 23, respectively. In addition, a 
comparison is made in figures 2k and 25 for a 12 -percent -thick Joukovski 
airfoil computed by the method of reference 3> and also by the method of 
this report. 

In examining figures 22 and 2k, which compare farthest position of 
impingement values, it can be seen that the calculated points compare 
fairly well with curves obtained as a result of formal solutions of the 
trajectory equations. A favorable comparison also appears to exist in 
the case of efficiency of impingement, as is shown by the curves in fig- 
ures 23 and 25. Discrepancies exist, but these may be expected in view 
of the empirical relations used as a basis for the computations. It 
will be noticed that values of farthest position of impingement were not 
obtained for values approaching the. limit position of drop impingement 
on upper and lower surfaces. This omission is due to the fact that the 
values of drop acceleration approach zero as the farthest position of 
impingement approaches the limit values on both upper and lower surfaces. 
As a result, scale moduli corresponding to these positions of impinge- 
ment are difficult to obtain accurately because a small error in drop 
acceleration makes a large error in scale modulus . For all practical 
purposes, though, inability to calculate points near the limit values 
makes little difference because the conditions in this range are seldom 
realized in flight. 
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Since the empirically calculated values for both farthest position 
and efficiency of impingement are in reasonable agreement with the 
theoretical values, it appears to be justifiable to presume that calcu- 
lations can be made with reasonable accuracy for other airfoils, at 
least in the same angle -of-attack range. Extension to other airfoils 
of the method using the hodograph technique tacitly includes the as- 
sumption that the vertical components of the tangential drop velocities 
bear the same relation to the maximum vertical velocities of the air as 
established in figure 10 for the five airfoils of this study. 
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TABLE I.- RESULTS FROM DIFFERENTIAL ANALYZER STUDIES OF WATER-DROP 
IMPINGEMENT ON A 15—IER CENT— THICK SYMMETRICAL JOUKDWSKT AIRFOIL 

[C2=0; a=0°] 


Figure 

number 


By 

y 


1(c) 

1(d) 

1(d) 

1(e) 

1(e) 

Ke) 

1(e) 

1(e) 

1(e) 

1(f) 

1(f) 

1(g) 

Kg) 

Kg) 

Kg) 

1(g) 
1(g) 
1(h) 
1(h) 
1 ( 1 ) 
1 ( 1 ) 
1( j) 
1(J) 

l(k) 

1(4) 

KD 

KD 

l(m) 

l( m ) 

l(m) 

l(m) 

l(m) 


32 

32 

4 

4 

16 
16 
16 
l6 
16 
16 
64 
64 
256 
256 
256 
256 
256 
■ 256 
8 
8 
32 
32 
128 
128 
512 
512 
2048 
2048 
64 
64 
64 
64 
64 


512 

512 

2048 

2043 

32 

32 

128 

128 

128 

128 

128 

128 

512 

512 

2043 

2048 

2048 

2048 

2048 

2048 

8 

8 

32 

32 

128 

128 

512 

512 

2048 

2048 

8 

8 

8 

8 

8 



0.074. 
-.074 
.074 
-.074 
.072 
-.072 
.073 
-.073 
.070 
.045 
.020 
-.020 
-.045 
-.070 
.0655 
-.0655 
.058 
.o4o 
.020 
-.020 
— . 040 
-.058 
.059 
— 059 
.056 
-.056 
.0485 
-.0485 
.038 
-.038 
.025 
-.025 
.0255 
.018 
.008 
_.008 
-.018 


Lower 1 

Upper 1 

Lawer 1 

Upper 1 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lower 1 

Upper 1 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lowery- 

Upper 1 

Lower 1 

Upper y- 

Lowery- 

Upper 1 

Lowery- 

Upper 1 

Lower 1 

Upper 1 

Lowery- 

Upper 1 

Upper 

Upper 

Lower 

Lower 


s/c 

V v 

V v 

0.265 

1.0 

0 

-.265 

1.0 

0 

.273 

.991 

.004 

-.273 

. 997 

-.004 

.262 

.997 

.013 

-.282 

.997 

-.013 

.273 

1.0 

.012 

-.273 

1.0 

-.012 

.244 

1.005 

.023 

.068 

.99 

.013 

.021 

.982 

.009 

-.021 

.982 

-.009 

-.068 

.99 

-.013 

-.244 

1.005 

-.023 

.225 

1.004 

.043 

-.225 

1.004 

-.043 

.188 

1.007 

.092 

.058 

.949 

.069 

.023 

• 931 

.029 

-.023 

.931 

-.029 

-.058 

.949 

-.069 

-.188 

1.007 

-.092 

.197 

.994 

. .078 

-.197 

• 994 

-.078 

.185 

.992 

.089 

-.185 

.992 

-.089 

.150 

.989 

.149 

-.150 

.9 89 

-.149 

.108 

.941 

.225 

-.108 

.941 

■ -.225 

.072 

.856 

.349 

-.072 

.856 

-.349 

.078 

.870 

.321 

.031 

.693 

.192 

.010 

.698 

.061 

-.010 

.698 

-.061 

-.031 

.693 

-.192 


L Denotes tangential trajectories. 
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TABLE I.- CONCLUDED 




% • 


Surface 

s/c 

%/v 

▼dA . 

l(m) 

l(n) 

l(n) 

l(o) 

l(o) 

l(o) 

l(o) 

l(o) 

l(o) 

l(p) 

1(P) 

l(q) 

l(q) 

Kr) 

l(r) 
l(s) . 
l(s) 
l(t). 
l(t) 

64 

256 

256 

1024 

1024 

1024 

1024 

1024 

1024 

4096 

4096 

16384 

16384 

512 

512 

8192 

8192 

32768 

32768 

8 

32 

32 

128 

128 

128 

128 

128 

128 

512 

512 

2048 

2048 

8 

8 

128 

128 

512 

512 

-0.0255 
.021 
-.021 
.015 
.010 
.005 
-.005 
-.010 
-.015 
.0110 
-.0110 
. 004 
-.004 

.0035 

-.0035 

.0020 

-.0020 

.0005 

-.0005 

Lower 1 

Upper 1 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower j 

Lower 1 

Upper 1 

Lawe^ 1 

Upper 1 

Lcwer 1 

Upper 1 

Lower 1 

Upper 1 

Lower 1 

Upper 1 

Lower 1 

-0.078 
.073' 
-.073 
.052 
.020 
.009 
-.009 
-.020 
-.052 | 

.038 
-.038 
.022 
-.022 
.023 
-.023 
.015 
-.015 
.016 
-.016 

: 0.870 
' .828 
.828 
.741 
,572 
.563 
.563 
.572 
.741 
,584 
.584 
.329 
.329 
.355 
.355 
.251 
.251 
.187 
.187 

-0.321 
.359 
-.359 • 
.451 
.198 
.109 
-.109 

-.198 
-. 451 
.452 
-.452 
.469 
-.469 
.514 
-. 514 
,4oi 
— 4oi 
.459 
-.459 

_ 


Denotes tangential trajectories. 
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TABLE II.- RESULTS FROM DIFFERENTIAL ANALYZER STUDIES OF WATER-DROP 
IMPINGEMENT ON A 15— PERCENT-THICK SYMMETRICAL JOUKDWSKI AIRFOIL 

[c z=0. 22, a=2° ] 


Figure 

number 


2(a) 

2(a) 

2(b) 

2(b) 

2(c) 

2(c) 

2(d) 

2(d) 

2(d) 

2(d) 

2(d) 

2(d) 

2(e) 

2(e) 

2(e) 

2(e) 

2(e) 

2(e) 

2(f) . 

2(f) 

2(f) 

2(f) 

2(f) 

2(f) 

2(g) 

2(g) 

2(h) 

2(h) 

2(h) 

2(h) 

2(h) 

2(h) 

2 (i ) 
2 ( 1 ) 
2(i) 
2(i) 


256 

256 

1024 

1024 

16 

16 

64 

64 

64 

64 

64 

64 

256 

256 

256 

256 

256 

256 

1024 

1024 

■1024 

1024 

1024 

1024 

16 

16 

64 

64 

64 

64 

64 

64 

256 

256 

256 

256 


-0.0046 

-.1548 

-.0055 

-.1539 

-.0081 

-.1533 

-.0095 

-.0381 

-.0667 

-.0956 

-.1243 

-.1532 

-.Ol4o 

-.o4io 

-.0683 

-.0958 

-.1232 

-.1508 

-.0214 

-.0464 

-.0721 

-.0977 

-.1231 

-.1488 

-.0435 

-.1598 

-.0493 

-.0705 

-.0902 

-.1130 

-.1345 

-.1558 

-.0587 

-.0763 

-.0940 
-. 1118 


Surface 


Upper 1 

Lower 1 

Upper 1 

Lowery- 

Upper 1 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lowery- 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lower 1 

Upper 1 

Lower 1 

Upper 1 

Upper 

Lower 

Lower 

Lower 

Lower 1 




Lower 


.223 
-.310 
. .226 
-.311 
.212 
■ .045 
.005 
-.026 
-.082 
-. 308 
.196 
.04l 
.003 
-.026 
-.078 
-.295 
.168 
.034 
.002 
-.027 
-.079 
-.265 
.149 
-.245 
.128 
.027 
-.001 
-.027 
-.071 
-.225 
.100 
.022 
-.002 
-.023 




1.001 

.998 

1.003 

.998 

1.009 

.997 

1.011 

.983 

.984 
' . 97 ^ 
.972 
.997 
1.015 
.980 
.969 
.970 
.975 
.995 
1.021 
.958 
.941 
.939 
.955 
.992 

1.010 

.984 

1.012 

.908 

.881 

.881 


0.04i 


^Denotes tangential trajectories. 
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21 
21 
21 
21 
21 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 

2(p) 

2(p) 

2(q) 

2(q) 



R v 

Jo 

Surface 

s/c 

%/7 

v d /v 

256 

256 

-O.1298 

Lower 

-0.058 

O.863 

-0.050 

256 

256 

-.1475 

Lower 1 

-.177 

.963 

-. 105 

1024 

1024 

-.0743 

Upper 1 

.063 

.914 

.424 

1024 

1024 

-.1360 

Lower 1 

-.118 

.902 

-.195 

128 

16 

-.0955 

Upper 1 

.055 

.881 

.472 

128 

16 

-.1065 

Upper 

.009 

.640 

.235 

128 

16 

-.1160 

Lower 

-.004 

.637 

.113 

128 

16 

-.1250 

Lower 

-.017 

.642 

-. 006 

128 

16 

-.1335 

Lower 

-.036 

.685 

-.109 

128 

16 

-.1430 

Lower 1 

-.098 

.854 

-.246 

512 

64 

-.1005 

Upper 1 

.038 

.767 

.528 

512 

64 

-.1080 

Upper 

.009 

.560 

.270 

512 

64 

-.1155 

Lower 

-.004 

.542 

.132 

512 

64 

-.1230 

Lower 

-.016 

.535 

-.011 

512 

64 

-.1310 

Lower 

-.03 

.611 

-. l4o 

512 

64 

-.1382 

Lower 1 

-.079 

.813 

-.300 

2048 

256 

^.1085 

Upper 1 

.028 

.562 

. .611 

2048 

256 

-.1130 

Upper 

.004 

.394 

.255 

2048 

256 

-.1182 

Lower 

-.005 

.379 

.118 

2048 

256 

1228 

Lower 

-.014 

.377 

-.050 

2048 

256 

-.1275 

Lower 

-.029 

.500 

-.235 

2048 

256 

-.1325 

Lower 1 . 

-.055 

.688 

-.356 

8192 

1024 

-.1165 

Upper 1 

.015 

.210 

.514 

8192 

1024 

-.1190 

— . 

0 

.270 

.230 

.8192 

1024 

1218 

Lower 

0 

0 

• 

1 

.260 

.085 

8192 

1024 

-.1250 

Lower ! 

-.012 

.273 

-. no 

8192 

1024 

-.1275 

Lower 1 

-.031 

.520 

-.390 

1024 

16 

-.1232 

Upper 1 

.008 

.186 

.610 

1024 

16 

-.1262 

Lower 1 

-.015 

.100 

-.243 

4096 

64 

-.1243 

Upper 1 

.003 

.090 

.434 

4096 

64 

-.1278 

Lower 1 

-.015 

.246 

-.235 

16384 

256 

-. 1254 ' 

Upper 1 

.002 

.103 

.506 

16384 

256 

-.1275 

Lower 1 

-.014 

.175 

.295 


■’’Denotes tangential trajectories. 
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TABLE III.- RESULTS FROM DIFFERENTIAL ANALYZER STUDIES OF WATER-DROP 
IMPINGEMENT ON A 1>-PERCENT-THICK SYMMETRICAL JOUKDWSKI AIRFOIL 

[c 2=0. 44; <x=4°] 


Figure 

number 

♦ 

% 


Surface 

s/c 

UdA 


EBB 

4 

256 

-0.1682 

Upper 1 

0.204 

0.9996 

O.O785 

SlOl. 

4 

256 

-.3215 

Lower 1 

-.4o8 

-.9946 

.0705 

3(4) 

l6 

1024 

-.1692 

Upper 1 

.194 

1.0086 

.1005 

3(4) 

16 

1024 

-.3223 

Lowery- 

-.400 

.9976 

.0705 

3(c) 

2 

16 

-.1818 

Upper 1 

.192 

1.0034 

.0962 

3(c) 

2 

16 

-.3330 

Lowery- 

-.409 

.9891 

.0688 

3(d) 

8 

64 

-.1837 

Upper 1 

.170 

1.0055 

.1032 

3(d) 

8 

64 

-.3073 

Lower 

-135 

.9802 

.0728 

3(d) 

8 

64 

-.2824 

Lower 

-.064 

.9782 

.0789 

3(d) 

8 

64 

-.2577 

Lower 

-.024 

.9663 

.1120 

3(d) 

8 

64 . 

-.2330 

Upper 

.004 

.9793 

.0940 

3(d) 

8 

64 

-.2083 

Upper 

.037 

.9854 

.1011 

3(d) 

8 

64 

-.3320 

Lowery- 

-.400 

.9930 

.O638 

3(e) 

32 

256 

-.1881 

Upper 1 

.148 

1.0174 

.1381 

3(e) 

32 

256 

-.2121 

Upper 

.034 

.9764 

.1221 

3(e) 

32 

256 

-.2358 

Upper 

.002 

.9653 

.1050 

3(e) 

32 

256 

-.259^ 

Lower 

-.024 

— 

— 

3(e) 

32 

256 

-.2832 

Lower 

-.062 

— 

— 

3(e) 

32 

256 

-.3068 

Lower 

, -.124 

.9652 

.0628 

3(e) 

32 

256 

-.3316 

Lower 1 

-.379 

.9831 

.0618 

3(f) 

128 

1024 

-.199* 

Upper 1 

.127 

1.0294 

.2031 

3(f) 

128 

1024 

-.2074 

Upper 

.062 

.9804 

.0931 

3(f) 

128 

1024 

-.2205 

Upper 

.028 

• 9514 

.1670 

3(f) 

128 

1024 

-.2386 

Upper- 

.004 

.9304 

.1410 

3(f) 

128 

1024 

-.2748 

Lower 

-.04l 

.9273 

.0869 

3(f) 

128 

1024 

-.3077 

Lower 

-.125 

• 9^33 

.0519 

3(f) 

128 

1024 

-.3316 

Lower 1 

-.350 

1.0200 

.0468 

3(g) 

16 

16 

-.2403 

Upper 1 

.121 

1.0116 

.2241 

3(g) 

16 

16 

-.3665 

Lower 1 

-.336 

.9628 

.0415 

3(4) 

64 

64 

-.2472 

' Upper 1 

.100 

1.0133 

.2881 

3(4) 

64 

64 

-.3422 

• Lower 

-.113 

.8940 

.0448 

3(4) 

64 

64 

-.3231 

Lower 

-.060 

.8712 

.0731 

3(4) 

64 

64 

-.3043 

Lower 

-.028 

.8694 

.1147 

3(4) 

64 

64 

-.2853 

Lower 

-.004 

.8638 

.1675 

3(4) 

64 

64 

-.2665 

Upper 

.018 

.8990 

.2148 

3(4) 

64 

64 

— .3606 

Lower 1 

-.286 

• 9578 

.0135 

3(i) 

256 

256 

-.2622 

Upper 1 

.068 


.4088 

3(i) 

256 

256 

-.2775 

Upper 

.012 

.8138 

.2746 

3(i) 

256 

256 

-.2925 

Lower 

-.008 

.7827 

.1843 


^■Denotes tangential trajectories. 













48 


NACA TN 2476 


TABLE III.- CONCLUDED 


Figure 

number 

* 

B V 

3b 

Surface 

. _ . 

s./c 

UdA 

mm 

wmwm 

-0.3078 

-.3353 

Lower • 
Lower 

-0.028 

-.084 

0.7864 

.8281 

3(i) 

256 

256 

-.3444 

Lower 

.-•125 

.8660 

3(1) 

256 

256 

-.3537 

Lower* 

-.247 

.9558 

3(j) 

1024 

1024 

-.2782 

Upper* 

.043 

.8758 

3(j) 

1024 

1024 

-.3440 

Lower* 

-.155 

.8780 

3(k) 

128 

16 

-.3126 

Upper* 

.042 

.8255 

3(k) 

128 

16 

-.3598 

Lower 

-.667 

• 71^3 

3(k) 

128 

16 

-.3500 

Lower 

— .040 

.6586 

3(k) 

128 

16 

-. 3406 

Lower 

— .022 

.6048 

3(k) 

128 

16 

-.3313 

Lower 

-.006 

.615 

3(k) 

128 

16 

-.3220 

Upper 

-.006 

.6562 

3(k) 

128 

16 

-.3688 

Lower* 

-.145 

.8632 

3(1) 

512 

64 

.-.3216 

Upper* 

.026 

.6392 

3(1) 

512 

64 

-.3303 

Lower 

-.004 

.5021 

3(1) 

512 

64 

-.3383' 

Lower 

-.015 

.4678 

3(1) 

512 

64 

-.3441 

Lower 

-.027 

.5197 

3(1) 

512 

64 

-.3500 

Lower 

-.042 

.5566 

3(1) 

512 

. 64 

-.3558 

Lower 

-.062 

.6465 

3(1) 

512 

64 

-.3605 

Lower* 

-.112 

. .8133 

3(m) 

2048 

256 

-.3293 

Upper* 

.015 

.4240 

3(m) 

2048 

256 

-.3324 

Upper 

.002 

.3300 

3 Cm) 

2048 

256 

-.3373 

Lower 

-.008 

.2858 

3(m) 

2048 

256 

-.3^32 

Lower 

-.022 

.35^8 

3(m) 

2048 

256 

-.3471 

Lower 

-.032 

.1+137 

3(m) 

2048 

256 

-.3501 

Lower • 

-.045 

.4956 

3(m) 

2048 

256 

-.3529 

Lower* 

-.072 

.6575 

3(a) 

8192 

1024 

-.3382 

Upper* 

.004 

.0918 

3(a)' 

8192 

1024 

-.3405 

Lower 

-.005 

.1918 

3(a) 

8192 

1024 

-.3^1 

Lower 

-.012 

.1787 

3(a) 

8192 

1024 

-.3458 

Lower 

-.020 

.2767 

3(a) 

8192 

1024 

-.3^ 

Lower 

-.031 

.2867 

3(a) 

8192 

1024 

-.3480 

Lower* 

-.041 

.4186 

3(o) 

1024 

16 

.-•3495 

Upper* 

.004 

.1147 

3(o) 

1024 

16 

-.35^4 

Lower 1 

-.035 

.31+96 

3(p) 

4096 

64 

-3^95 

Upper* 

.002 

.1047 

3(p) 

4096 

64 

-.3530 

Lower* 

-.026 

.2146 

3(q) 

16384 

256 

-.3504 

Upper* 

: .002 

.1597 

3(q) 

16384 

256 

-•3535 

Lower* 

-.025 

.2746 


v d /v 

0.1092 
.0149 
.0142 
^.0313 
■ .5817 
-.1212 
• 5550 
.0364 
.0228 
.1031 
.2444 
. 3^37 
-.1358 
.7017 
.3024 
.1513 
.0390 
-.0782 
-.1363 
-.2045 
.7513 
A553 
.2442 
.0060 
-.1191 
-.1923 
-.2912 
.7521 
.4371 

.1520 

.0669 

-.2011 

-.2832 

.6545 

-.2808 

.6945 

-.2157 

.8494 


denotes tangential trajectories. 
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TABLE IV.- RESULTS FROM DIFFERENTIAL ANALYZER STUDIES OF WATER-DROP 
IMPINGEMENT ON A 15 — PERCENT-THICK CAMBERED JOUKOWSKI AIRFOIL 
[a= 1.0 MEAN LINE; c 2 = 0 . 44 ; a= 0 ° ] 


Figure 

number 

* 

% 


Surface 

s/c 

m/v 

V 7 - 

WB 3 M 

4 

256 

0.0935 

Upper 1 

0.317 

1.003 

0.007 

■sm 

4 

256 

-.0565 

Lover 1 

-.216 

.998 

-.002 

E wm 

16 

1024 

.0915 

Upper 1 

.310 

1.008 

.013 

Mb) 

16 

1024 

-.0565 

Lover 1 

-.213 

• 999 

-.006 

Me) 

2 

16 

.0855 

Upper 1 

.305 

1.009 

.022 

Me) 

2 . 

16 

-.0600 

Lover 1 

-.215 

.994 

-.001 

Md) 

8 

64 

.0818 

Upper 1 

.294 

. 1.012 

.031 

Md) 

8 

64 

.0536 

Upper 

.096 

.992 

.022 

Md) 

8 

64 

.0253 

Upper 

.038 

.988 

.020 

Md) 

8 

64 

-.0038 

Upper 

.002 

.985 

.011 

Md) 

8 

64 

-.0318 

Lower 

-.033 

.984 

.003 

Md) I 

8 

64 

-.0610 

Lower 1 

-.212 

.999 

-.008 

Me) 

32 

256 

.0775 

Upper 1 

.275 

1.022 

.050 

Me) 

32 

256 

.0503 

Upper 

.092 

.989 

.042 

4 (e) 

32 

256 

.0225 

Upper 

.034 

.976 

.030 

4 (e) 

32 . 

256 

-.0045 

Upper 

.001 

.972 

.017 

Me) 

32 

256 

-.0325 

Lower 

-.033 

.973 • 

.002 

Me) 

32 

256 

-.0600 

Lover 1 

-.199 

.990 

-.016 

4 (f) 

128 

1024 

.0660 

Upper 1 

.243 

1,033 

.090 

4 (f) 

128 

1024 

.0420 

Upper 

.077 

.976 

.075 

Mf) 

128 

1024 

.0160 

Upper 

.030 

•954 

.054 

Mf) 

128 

1024 

-.0085 

— 

0 

.943 

.024 

Mf) 

128 

1024 

-.0335 

Lover 

-.032 

.946 

-.006 

4 (f) 

128 

1024 

-.0585 

Lover 1 

-.183 

.984 

-.036 

Mg) 

16 

16 

.0377 

Upper 1 

.211 

1.028 

.128 

Mg) 

16 

16 

-.0770 

Lover 1 

-.180 

.978 

-.042 

Mh) 

64 

64 

.0312 

Upper 1 

.192 

1.038 

.168 

4 (h) 

64 

64 

.0100 

Upper 

.061 

.936 

.131 

4 (h) 

64 

64 ' 

-.0110 

Upper 

.022 

.898 

• 095 

Mh) 

64 

64 

-.0315 

Lover 

-.003 

.884 

.048 

Mh) 

64 

64 

-.0525 

Lover 

-.031 

.891 

-.007 

4 (h) 

64 

64 

-.0731 

Lover 1 

-.157 

.971 

-.073 

Mi) 

256 

256 

.0180 

Upper 1 

.158 

1.038 

.238 

Mi)’ 

256 

256 

.0010 

Upper 

.050 

• 893 

.188 

Mi) 

256 

256 

-.0165 

Upper 

.018 

.839 

.127 

Mi) 

256 

256 

-.0340 

Lover 

-.004 

.820 

.054 

MD 

256 

256 

-.0510 

Lover 

-.028 

.838 

-.017 

Mi) 

256 

256 

-.0680 

Lover 

-.120 

.940 

-.129 


L Denotes tangential trajectories. 
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TABLE IV.- CONCLUDED 


Figure 

number 

If 

% 

y 0 

Surface 

s/c 

Ud/V 


Mj) 
Mj) 
Mk) 
Mk) 
Mk) 
Mk) 
Mk) 
4(k) 
I Ml) 
Ml) 
MD 
MD 
MD 
MD 

Mm) 

Mm) 

Mm) 

Mm) 

Mm) 

4(m) 

Mn) 

Mn) 

Mn) 

Mn) 

Mn) 

Mn) 

Mo) 

Mo) 

Mp) 

Mp) 

Mq) 

Mq) 

1024 

1024 

128 

128 

128 

128 

128 

128 

512 

512 

512 

512 

512 

512 

2048 

2048 

2048 

2048 

2048 

2048 

8192 

8192 

8192 

8192 

8192 

8192 

1024 

1024 

4096 

4096 

16384 

16384 

'1024 

1024 

16 

16 

16 

16 

16 

16 

64 

64 

64 

64 

64 

64 

256 

256 

256 

256 

256 

256 

1024 

1024 

1024 

1024 

1024 

1024 

16 

16 

64 

64 

256 

256 

0 

-.0620 

-.028 

-.0382 

-.0480 

-.0575 

-.0670 

-.0772 

-.0347 

-.0425 

-.050 

-.0582 

—.066 

-.0725 

-.0440 

-.050 

-.0548 

-.0594 

—.0638 

-.0685 

-.0548 

-.0569 

-.0590 

-.0610 

-.0630 

-.0650 

-.0604 

-.0700 

-.0650 

-.0675 

-.0655 

-.0670 

Upper 1 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lower 1 

Upper 1 

Upper 

Upper 

Lower 

Lower 

Lower 1 

Upper 1 

Upper 

Lower 

Lower 

Lower 

Lower 1 

Upper 1 

Lower 1 

Upper 1 

Lower 1 

Upper 1 

Lower 1 

0.109 

-.085 

.092 

.032 

.012 

-.004 

-.019 

—.068 

.072 

.025 

.010 

-.005 

-.018 

-.053 

.046 

.015 

.005 

-.004 

-.013 

-.038 

.025 

.007 

-.001 

-.008 

-.012 

-.025 

.022 

-.018 

.015 

-.012 

.006 

-.008 

1.000 

.883 

.958 
.727 
.6 38 
.649 . 
.644 
.838 
• 911 
.634 
.535 
.485 
.544 
.754 
.686 
.448 
.381 
.368 
.407 
.604 
.471 
.258 
.178 
.184 
.211 
.464 
.346 
.317 
.285 

.143 

.113 

.046 

0.357 
-.220 
.405 
.288 
.189 
.060 
-.046 
-.264 
.487 
. -350 
.208 
.031 
-.124 
-.319 
• 576 

.360 

.192 

.048 

-.103 

-.384 

.650 

.308 

.082 

-.079 

-.168 

-.428 

.561 

-.368 

.527 

-.253. 

.391 • 
-.118 


^De notes tangential trajectories. 
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TABLE V.- RESULTS FROM DIFFERENTIAL ANALYZER STUDIES OF 
WATER-DROP IMPINGEMENT ON AN NACA 65 2 -015 AIRFOIL 

[0^=0. 44; a=4°] 


Figure 

number 

♦ 

m 


Surface 

s/c 

udA 

v d A 

5(a) 

4 


-0.1281 

Upper 1 

0.281 


0.0814 

5(a) 

4 


-.2817 

Lover 1 

-.523 


.0704 

5(b) 

16 


-.1298 

upper 1 

.267 

1.0073 

.0864 

5(b) 

16 

1024 

-.2818 

Lower 1 

-.514 

.9973 

.0704 

5(c) 

2 

16 

-.1395 

Upper 1 

.259 

1.0047 

.0945 

5(c) 

2 

16 

-.1646 

Upper 

.031 

.9881 

.0875 

5(c) 

.2 

16 

-.2026 

Lower 

-.022 

.9847 

.0825 

5(c) 

2 

16 

-.2246 

Lower 

-.074 

.9867 

.0793 

5(c) 

2 

16 

-.2467 

Lower 

-.150 

.9857 

•0753 

5(c) 

2 

16 

-.2687 

Lower 

-.256 

.9847 

.0742 

5(c) 

2 

16 

-.2909 

Lower 1 

-.514 

.9927 

.0711 

5(d) 

8 

64 

-.1424 

Upper 1 

.240 

1.0107 

.1065 

5(d) . 

8 

64 

-.1719 

Upper 

.016 

.9847 

.0935 

5(d) 

8 

64 

-.2163 

Lower 

' -.050 

.9757 

.0824 

5(d) 

8 

64 

-.2409 

Lower 

-.125 

.9807 

.0743 

5(d) 

- 8 

64 

-.2655 

Lower 

--.236 

.9797 

■ .0712 

5(d) 

8 

64 

-.2899 

Lower 1 

-.512 

.9907 

.0641 

5(e) 

32 

256 

-.1493 

Upper 1 

.209 

1.0187 

.1326 

5(e) 

32 

256 

-.1702 

Upper 

.023 

.9761 

.1155 

5(e) 

32 

256 

-.2193 

Lower 

-.052 

.9617 

.0874 

5(e) 

32 

256 

-.2437 

Lower 

-.131 

• 9647 

.0743 

5(e) 

32 

256 

-.2685 

Lower 

-.249 

.9717 

.0632 

5(e) 

32 

256 

-.2894 

Lower 1 

-.506 

1.0057 

.0531 

5(f) 

128 

1024 

-.1603. 

Upper 1 

.150 

1.0207 

.1805 

5(f) 

128 

1024 . 

-.1826 

Upper 

.008 

.9517 

.1375 

5(f) 

128 

1024 

-.2282 

Lower 

-.068 

.9427 

.0853 

5(f) 

128 

1024 

-.2505 

Lower 

-.145 

.9407 

.0642 

5(f) 

128 

1024 

-.2726 

Lower 

-.267 

• 9557 

.0462 

5(f) 

128 

1024 

-.2878 

Lower 1 

-.485 

.9837 

.0382 

5(g) 

16 

16 

-.1951 

Uppen 

.128 

1.0015 

.2008 

5(g) 

16 

16 

-.3202 

Lowen 

-.481 

.9776 

.0383 

5(b) 

6k 

64 

-.2030 

Upperi 

.092 

.9956 

.2597 

5(b) 

6k 

64 

-.2136 

Upper 

.012 

.9155 

.2055 

5(b) 

6k 

64 

-.2535 

Lower 

-.048 

.8795 

.1092 

5(b) 

6k 

64 

-.2787 

Lower 

-.127 

.8975 

.0609 

5(h) 

6k 

64 

-.3013 

Lower 

-.247 

.9255 

.0227 

5(h) 

6k 

64 

-.3111 

Lowery- 

-.417 

.9756 

.0074 

5(i) 

256 

256 

-.2143 

Upper 1 

.047 

.9416 

.3416 

5(1) 

256 

256 

-.2343 

Lower 

-.007 

.8256 

.1964 


1 Denotes tangential trajectories. 
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TABLE V.- CONTINUED 


Figure 

number 

♦ 

By 

7o 

Surface 

s/c 


VdA 

5(i) 

256 

256 

-0.2543 

Lower 

-0.040 

0 . 8086 

0.1102 

5(i) 

256 

256 

-.2717 

Lower 

-.093 

.8316 

.0410 

5(i) 

256 

256 

-.2892 

Lower 

-.184 

— 

— 

5(i) 

256 

256 

-.2983 

Lower 1 

-.355 

.9686 

-.0443 

5(j) 

1024 

1024 

-.2267 

Upper 1 

.026 

.8575 

.4614 

5( j) 

1024 

1024 

-.2313 

Upper 

.010 

.7726 

.3935 

5(j) 

1024 

1024 

-.2492 

Lower 

-.015 

.6715 

.1972 

5(j) 

1024 

1024 

-.2692 

Lower 

-.060 

.6295 

.0300 

5( j) 

1024 

1024 

-.2782 

Lower 

-.104 

.7875 

-.0251 

5( j) 

1024 

1024 

-.2853 

Lower 

-.166 ' 

.7^95 

-.0692 

5(j) 

1024 

1024 

-.2883 

Lower 1 

-.245 

.9175 

.0942 

5(k) 

128 

16 

-.2633 

Upper 1 

.021 

.7926 

.4950 

500 

128 

16 

-.2798 

Lower 

-.012 

.6436 

.2336 

5(k) 

128 

16 

-.2945 

Lower 

-.042 

.6395 ' 

.1003 

5(k) 

128 

16 

-.3028 

Lower 

-.080 

.7194 

-.0500 

5(k) 

128 

16 

-.3059 

Lower 

-.104 

.7634 

-.0860 

5(k) 

128 

16 

-.3091 

Lower 1 

-.185 

.8804 

-.1302 

5(1) 

512 

64 

-.2676 

Upper 1 

.018 

.7536 

.5779 

5(1) 

512 

64 

-.2719 

Upper 

.005 

.6496 

.4378 

5(1) 

512 

64 

-.2857 

Lower 

-.020 

.^935 

.1605 

5(1) 

512 

64 

-.2926 

Lower 

-.034 

.5504 

.040^ 

5(1) 

512 

64 

-.3015 

Lower 

-.079 

.6954 

-.0989 

5(1) 

512 

64 

-.3035 

Lower 1 

-.145 

.8284 

-.1540 

5(m) 

2048 

256 

-.2737 

Upperi 

.014 

.6606 

.6718 

5(m) 

2048 

256 

-.2748 

Upper 

.007 

.5986 

. 5868 

5(m) 

2048 

,256 

-.2781 

— 

0 

.5216 

.4106 

5(m) 

2048 

256 

-.2861 

Lower 

-.012 

.4085 

.2245 

5(m) 

2048 

256 

-.2906 

Lower 

-.022 

.3905 

.0604 

5(m) 

2048 

256 

-.2965 

Lower 

-.048 

.5485 

-.0967 

5(m) 

2048 

256 

-.2989 

Lowery- 

-.100 

.7584 

-.1768 

5(n) 

8192 

1024 

-.2798 

Upper 1 

.010 

.5976 

.7847 

5(n) 

8192 

1024 

-.2801 

Upper 

.005 

— 

— 

5U) 

8192 

1024 

-.2850 

Lower 

-.004 

— 

— 

5(n) 

8192 

1024 

-.2886 

Lower 

-.010 

.3035 

.2454 

5(n) 

8192 

1024 

-.2932 

Lower 

-.022 

.2124 

.0033 

5(n) 

8192 

1024 

-.2945 

Lower 

-.030 

.3423 

-.0868 

5(n) 

8192 

1024 

-.2958 

Lower 

-.033 

.3983 

-.1068 

5(n) 

8192 

1024 

-.2971 

Lowery- 

-.065 

.6765 

-.1867 

5(o) 

1024 

16 

-.2933 

Upper 1 

.008 

.4372 

.8238 

5(o) 

1024 

16 

-.3037 

Lower 1 

-.040 

.4891 

-.1997 


x Denotes tangential trajectories. 
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TABLE V.- CONCLUDED 
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TABLE VT.- VALUES OF C d R/24 AS A FUNCTION OF R 


- - R 

C d R/24 

R 

C d R/24- 

0.00 

1.000 

500 

11.46 

• .05 

1.009 

600 

12.97 

.1 

1.018 

800 

15.81 

- .2 

1.037 

1,000 

18.62 


1.073 

1,200 

21.3 

.6 

1.108 

1,400 

24.0 

.8 

1.142 

1,600 

26.9 

- 1.0 

1.176 

1,800 

29.8 

-1.2 

1.201 

2,000 

32.7 

1.4 

1.225 

2 -,500 

40.4 

1.6 

1.248 

3,000 

47.8 

1.8 

1.267 

3,500 

55.6 

■2.0 

1.285 

4,000 

63.7 

2.5 

1.332 

5,000 

80.0 

3.0 

1.374 

6,000 

'96.8 

3.5 

1.412 

8,000 

130.6 

4.0 

1.447 

10,000 

166.3 

5.0 

1.513 ; 

12,000 

204 

.6.0 

■ 1.572 

14,000 

243 

8.0 

■ 1.678 

16,000 

285 

10.0 

1.782 

18,000 

325 

12 

1.901 

20,000 

365 

14 

2.008 

25,000 

470 

16 

2.109 

30,000 

574 

18 

2.198 

35,000 

674 

20 

2.291 

40,000 

778 

25 

2.489 

50,000 

980 

30 

2.673 

60,000 

1175 

35 

2.851 

80,000 

1552 

- 4o 

3.013 . 

100,000 

1905 

' 50 

3.327 ' 

1.2x105 

2234 

- 60 

3.60 

1.4X105 

2549 

■ 80 

4.11 

1.6X10 5 

2851 

100 

4.59 



120 

5.01 ■ 



140 

5.40 



■ 160 

5.76 



180 

6.16 



20 b 

1 6.52 



250 

7.38 ■ 



300 

8 .26 



350 

9 .00 



400 

9.82 








TABLE VII.- COMPARISON OF PRODUCTS CF SCALE MODULUS AND DROP DRAG COEFFICIENT FOR THE 
UPPER AND LOWER SURFACES OF A 15-PERCENT-THICK SYMMETRICAL JOUKDWSKI AIRFOIL 

[a=0°j cj=0] 


NACA TN 2476 


PvTI 
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ooh ro 
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ON ON-4 
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H 
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01 
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AN NACA 65 2 -015 AIRFOIL 
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(a) !5 - percent- thick symmetrical Joukowski ; 




Upper surface 



(b) 15- percent -thick symmetrical Joukowski ; 





Upper surface 
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Figure /. - Concluded. 
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Figure 2- Concluded. 
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Figure 6-Water-drop trajectory traces from a differential analyzer for an NACA 65-0/5 
airfoil ■, » 0.44 1 a = 4°. ‘ 


y, percent y , percent \, \ | | y , percent \ ( \ \ y t percent 
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Figure 6 -Continued. 
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Figure 6. - Continued. 
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Scale modulus, 

Figure 7 - Typical relation between farthest position of drop impingement, scale modulus, and 
free -stream drop Reynolds number as exemplified by a 15- percent -thick cambered 

Joukowski airfoil; c= 0.44; a=0° ; a = 1.0 mean line. 
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u„ u d 

x- velocity component of air ; —, and of drop, — 


(a) 15 -percent -thick symmetrical Joukowski airfoil, a - 0°. 

Figure 8.- Hodograph of tangential -trajectory velocities and 
hodograph of air velocities on the airfoil surface. 
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(b) 15 - percent -thick symmetrical Joukowski airfoil a = 2°. 
Figure 8. —Continued . 
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x -velocity component of air, -y-, and of drop, 


(c) !5- percent -thick symmetrical Joukowski airfoil, a *4° 


Figure 8. — Continued. 
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(d) 15 -percent -thick cambered Joukowski airfoil; a = 1.0, 
a -0°. 


Figure 8 . - Continued. 
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(e) NACA 65 r 0!5 airfoil, a ‘4°. 


Figure 8..- Concluded . 
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(a) Physical plane . 



(b) Hodograph plane . 



Figure 9. - Relationship between physical and hodograph planes for 
drop and air velocities at airfoil surface . 
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Figure 10.— Variation of velocity difference between drop and air with maximum 
y -velocity component of air for the five airfoil cases investigated . 



Curve C based on » constant 
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x- velocity component of air, y , and of drop, -y 

Figure //.- Illustration of four possible techniques for the construction of a drop 
hodograph from a specified air hodograph. 
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Drop acceleration 






Drop acceleration 
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(d) 15- percent - thick cambered Joukowski 
airfoil ; a - 0°; a - i.O mean line . 

Figure 12,- Continued . 





Drop acceleration 


NACA TN 2476 


105 




0 10 20 30 40 50 

Chord wise position , x, percent 

(e) NACA 65 g -0/5 airfoil ; a - 4° . 


Figure !2. - Concluded . 
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(b) 45- percent- thick symmetrical Joukowski airfoil ; 
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Figure ! 3 - Concluded. 
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2 .3 .4 .5 

Lift coefficient , c, 


(a) Joukowski airfoil. 


Figure 14. Velocity gradient along the stagnation streamline, at 
the stagnation point, as a function of lift coefficient and 


thickness ratio for two profiles. 
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(b) Ellipse. 


Figure 14. Concluded. 







Ratio of starting ordinate Airfoil ordinate at farthest Distance along airfoil surface 
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Scale modulus, \jr 

(c) Ordinate- intercept ratio 

- Graphical representation of the procedure used to 
a value of \fr corresponding to( Ay ^ Ay ) t =0.8 
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(a) 15- percent -thick symmetrical Joukowski airfoil ; a - 0°. 

Figure 17 - Surface distribution of water-drop impingement for five 
airfoils. 
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(c) 15 -percent -thick symmetrical Joukowski airfoil ; 
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o -20 -to O 10 20 30 

Distance along airfoil surface from leading edge, s/c, percent 
(d) 15- percent -thick cambered Joukowski airfoil ; a - i.O mean line ; 
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Figure 17. — Concluded . 
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Figure 18 — Concluded. 




NAG A 23015 airfoil-, a .* 3.6 ° ; c.= 0.5 . 


Distance along airfoil surface from leading edge , s/c, percent 
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Figure 19. - Continued . 
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x, percent chord 

(c) Slope of airfoil contour as a function of chord wise position . 
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L L 

x- velocity component of air, -y- 
(e) Hodograph for air. 


Figure 19. -Continued. 
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(f) Initial phase of drop hodograph . 


Figure 19. -Continued. 




- velocity component of air, -j y , and of drop. 
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x -velocity component of air, — , and of drop, — 


(g) Completed drop and air hodographs . 
Figure 19. - Continued . 
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Calculated points 
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(0 Distribution of drop acceleration values over airfoil surface. 




Actual drag coefficient as a multiple of drag coefficient 
given by Stokes law, C d R/ 24, dimensionless 
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Figure 19. - Continued . 
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(k) Farthest position of impingement . 
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Figure 20. - Graphical relationships used in evaluating impingement efficiency 
for an NACA 230/5 airfoil, a =3.6°; c t -0.5 . 
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(c) Ordinate -ratio isop/eths. 
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(c) Distribution of impingement. 
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(a) 15-percent-thick symmetrical Joukowski airfoil ; a s 0°. 

Figure 22- Empirically calculated point values of farthest position of impingement for five 
airfoils in comparison with curves obtained from a differential analyzer. 
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(d) 15 -percent -thick cambered Joukowski airfoil ; a b 1.0 mean line; 
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Figure 22 - Concluded. 
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Figure 23. Empirically calculated point values of efficiency of impingement it 
comparison with curves obtained from a differential analyzer. 
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(b) 15-percent- thick symmetrical Joukowski airfoil ; a=2°. 
Figure 23.— Continued. 
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Figure 23 — Continued . 
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(d) 15- percent-thick cambered Joukowski airfoil; a *1.0 mean line -, 
Figure 23— Continued. 
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Figure 23 — Concluded. 






150 


NACA TN 2476 



fU90J9d ‘D/S *96p9 
6UfpD9j UJOJJ 9DDJJHS l/OJJID 5U0/0 9DUDJSIQ 


I 

I 

f 

V 

1 

<0 

0) 

$ 

I 

V . 

0 

5 

1 

1 

I 

t 

I 

I 

t 

I 


I 

1 


1 

* 

SM 


I 

I 


! 

t 

t 

<0 

•o' 

* 

I Q 

s - 

■§ * 

$ ¥ 
5 *5 

II 

I ^ 

$ t 
1 3 










NACA TN 2476 


151 



tudojsd *3 V uduidBufdujf jo Aoudjofjj^ 


% ! 

*1 

s I 

* * 

§ £ 

52 £ 

I 5 

1 1 . 

§ «■ 

s l 

> < 

I i 

II 

.1 ^ 

$ cj 

g* 

s £ 

0 .1 


s 

s 


i 


i 

I 

ji^ 


1 * 

: l 1 • 

I I ^ 

Uj I 1 
10 

<M <g 


.M 


5» < 


I 


NACA-Langley - 9-14-51 - 1000 


